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I. INTRODUCTION 

The statistical mechanics of hquids and fluids is by now a mature field which has evolved 
from the study of simple liquids to that of increasingly complex systems and phenomena 
(see e.g [1]). A research direction which has received only limited attention so far is that of 
liquids embedded in curved spaces. Although the topic might appear at first sight as a purely 
theoretical curiosity or only relevant to general relativity, it arises in situations in which 
particles are adsorbed or confined on a substrate with nonzero curvature, be it the wall of a 
porous material, the surface of a large solid particle or an interface in an oil-water emulsion. 
In the latter cases, the fluid is confined to a curved manifold which is two-dimensional. 
Statistical-mechanical systems in three-dimensional curved spaces and manifolds are not 
encountered in our everyday life experience, but they are useful as templates or models 
to study properties whose direct investigation in Euclidean three-dimensional space remains 
difficult or inconclusive: curvature and geometry then provide additional control parameters 
to envisage the behavior of a system. 

In this article, we review the progress made on the statistical mechanics of liquids and 
fluids embedded in curved space. Our main focus will be on two-dimensional manifolds 
of constant nonzero curvature and on the influence of the latter on the phase behavior, 
thermodynamics and structure of simple liquids. Reference will also be made to existing work 
on three-dimensional curved space and two-dimensional manifolds with varying curvature. 
On the other hand, we exclude from the scope of the article substrates with fluctuating 
geometry, such as membranes [21 [3]. The geometry will always be considered as frozen, 
providing the background whose metrical and topological characteristics affect the behavior 
of the embedded fluids, but with no feed-back influence from the latter. 

The rest of the paper is organized as follows: in Sec. [IT], we review physical examples where 
the curvature of substrate modify the physical properties of systems in comparison with those 



observed in in Euclidean space. In Sec. Ill, we discuss the speciflcities of thermodynamics for 



liquids coming from the flniteness of space (constant positive curvature) or strong influence 



of the boundary effects (constant negative curvature). We introduce in Sec. IV elements of 



liquid state theory for describing the structure and the thermodynamic properties of fluids 



in curves spaces. Section W\ is devoted to the influence of tlie curvature on tlie structure of 



liquids as well as the modifications of the Coulomb interaction We consider in Sec. VI ill- 
ordered dense phases where the geometric frustration prevents the appearance of an ordering 
transition and slows down drastically the relaxation dynamics of liquids driving the system 



to a glassy behavior. In Sec. VII, the low-temperature regions of phase diagram is analyzed 



by means the elastic theory of defects, which reveals the structure of topological defects 



generally in excess in curved spaces. Concluding remarks are drawn in Sec. VIII 



II. THEORETICAL MOTIVATIONS AND PHYSICAL REALIZATIONS 

In this section we expand a little more on the reasons for studying fluids in curved spaces. 
Before moving on to more solid grounds, let us first acknowledge that non-Euclidean geome- 
tries are fun and fascinating. Musing about the differences between spherical and hyperbolic 
worlds, i.e. behavior in spaces of constant positive and negative curvature respectively, is 
an exciting intellectual experience! At a fundamental level, a nonzero curvature introduces 
(at least) one extrinsic lengthscale (a "radius of curvature") in the behavior of systems em- 
bedded in such curved spaces. One then expects that the "long-distance" properties of the 
system, namely those involving lengths much larger than the radius (or radii) of curvature, 
are modified whereas "local" ones should be rather insensitive to curvature in general. Curv- 
ing space may also change the topology of the substrate, going for instance from an infinite 
flat plane to a sphere or a torus, which also affects some properties of the embedded system 
such as the nature of its ordered condensed phases. 

From a theoretical point of view, curving space provides an additional control parameter 
(the Gaussian curvature or the associated radius of curvature for a homogeneous space 
of constant curvature) for studying the properties of a fluid or a liquid, in addition to the 
common thermodynamic parameters. This may prove interesting in several situations. First, 
there are cases for which in the standard Euclidean space, boundary conditions matter. This 
is true for instance for Coulombic systems in which charged particles interact through a long- 
range Coulomb potential or in the vicinity of a gas-liquid critical point where correlations 
extend over the whole system size. The long-range character of the interactions or of the 
correlations entail the use of boundary conditions, usually periodic boundary conditions, 
in standard statistical-mechanical treatments. An alternative is provided by the so-called 



"spherical and hyperspherical boundary conditions" jlHB] , which amount to considering the 
system on the surface of a sphere (for a two-dimensional space) or a hypersphere (for a three- 
dimensional space). A (hyper)sphere is a homogeneous and isotropic manifold of constant 
positive Gaussian curvature. It is of finite extent so that no boundary conditions need to be 
specified. The Euclidean space is then recovered by letting the radius of the (hyper) sphere 
go to infinity. 

Another domain in which curving space proves to be of major theoretical importance is 
what goes under the concept of "geometric frustration" [7] . The latter describes an incom- 
patibility between the preferred local order in a system and the tiling of the whole space. 
Geometric frustration has emerged in the theoretical description of glasses and amorphous 
solids from the consideration of local icosahedral order in metallic glasses. Icosahedral, or 
more generally poly-tetrahedral order, is favored for local arrangement of atoms [8J but can- 
not extend to form a periodic tiling of three-dimensional Euclidean space [H [H]. It has first 
been realized by Kleman and Sadoc[Tn], and further developed by several groups[71 19], that 
perfect tetrahedral/icosahedral order can exist on the surface of a hypersphere with an ap- 
propriately chosen radius. More generally, curving space can be seen as a way to reduce (as 
in the above example) or increase the amount of geometric frustration in a liquid. Increasing 
frustration is helpful for liquids of spherical particles in two dimensions. Indeed, the ordinary 
"flat " space then leads to no frustration. The locally preferred hexagonal order can tile space 
to produce a triangular lattice. As a result, one-component atomic liquids on a Euclidean 
plane rapidly order under cooling (or compressing) and do not form glasses. To study glass 
formation and amorphous packings in such systems, one must introduce curvature. As sug- 
gested by Nelson and coworkers [9| [TT| [T2] . one could then mimic frustrated icosahedral order 
in ordinary three-dimensional space by considering frustrated hexagonal order in the hyper- 
bolic plane, which is a homogeneous and isotropic two-dimensional manifold of constant 
negative Gaussian curvature. 

Up to now, we have focused on purely theoretical motivations for studying curved spaces. 
As mentioned in the introduction, there are also physical realizations of fluids on two- 
dimensional curved substrates. Such situations occur in adsorption and coating phenomena. 
Generally, the adsorbing solid substrate is not fiat (it could be locally cylindrical or spheri- 
cal) and an equilibrium fluid monolayer can form on its surface if the fluid-solid attraction is 
strong enough. There are cases where the curvature can even be quite strong, with the asso- 



dated length only moderately larger than that of the adsorbate. This happens for instance 
in the adsorption of gas molecules inside the cavities of zeolite molecular sieves or that of 
colloidal particles on larger spherical particles in the phenomenon of heterofioculation|13j. 

Fluid and condensed-phase behavior is also observed at the oil- water interface in emulsions [TU 
[T3] . The presence of colloidal particles at the interfaces actually stabilize the emulsions, 
then often referred to as Pickering emulsions |16j. The colloidal particles are irreversibly 
adsorbed but still mobile at the interface which is usually the surface of a spherical droplet. 
Up to now, most of the existing experiments focus on dense, crystal-like arrangements 
of colloidal particles on spherical droplets [HI [151 [T7H20] . but one could imagine studying 
less dense, liquid-like behavior [21J. One could also consider adsorbed colloidal particles on 
more exotic interfaces such as the infinite periodic minimal surfaces formed for instance by 
amphiphilic molecules in the presence of water [22142^ . Such surfaces have negative, but 
varying, curvature and are periodic in all three directions. On the other hand, we do not 
discuss in this review situations in which the geometry of the substrate changes when coated 
by colloidal particles, as may occur for bubbles in foams for instance [25l [26] . In addition, 
we only briefly allude to crystalline-like order on curved substrates, as this important topic 
has been recently extensively reviewed by Bowick and Giomi[27]. 

III. THERMODYNAMICS AND BOUNDARY EFFECTS 

A. Thermodynamic limit and boundary conditions 

Thermodynamics, and more generally all dynamic and structural properties of a macro- 
scopic system, are retrieved with statistical mechanics by taking the thermodynamic limit 
in which one let the number of particles and the volume of the system go to infinity while 
keeping the density finite. This also ensures equivalence between the various thermody- 
namic statistical ensembles. Such a procedure is, at least conceptually, easy to implement in 
Euclidean space when both interactions and spatial correlations between particles are short- 
ranged. So long as the system is finite, the properties under study depend on the choice of 
the boundary conditions, e.g. free or periodic, but as one takes the thermodynamic limit 
the contribution of the boundary vanishes irrespective of the chosen conditions (the rate 
of convergence does however depend on the boundary conditions). The same behavior also 



characterizes manifolds of varying curvature that are of infinite extension and can be em- 
bedded in Euchdean space. Note however, that in addition to the subtleties occurring for 
long-range, e.g. Coulombic, interactions or long-range correlations near a critical point, some 
care must be exerted for anisotropic and multi-connected manifolds as the spatial extent of 
the system can be made infinite only in certain directions. 

Quite different is the behavior of homogeneous spaces of constant nonzero curvature. 
Spherical geometry associated with a positive curvature indeed leads to spaces of finite 
extent whereas hyperbolic geometry associated with a negative curvature allows for infinite 
spatial extent but with a strong effect of the boundary conditions. In what follows, we 
discuss these two main geometries, specializing for simplicity to the two-dimensional case. 

B. Spherical geometry 

In two dimensions, the space of constant positive curvature is the surface 5*2 of a sphere. 
If R is the radius of the sphere, l/R"^ is its Gaussian curvature and its area AnR"^ is obviously 
finite. Taking the thermodynamic limit i? — )■ oo "flattens" the space and one then recovers 
the Euclidean plane. The same is true in higher dimensions. If one wishes to consider fluids 
on spherical substrates of constant curvature l/i?^, one must resort to the thermodynamics 
of finite systems. This topic has been (ands still is) extensively discussed. For small systems, 
one route is provided by Hill's formulation [28] of the thermodynamics of small systems: one 
starts with a finite canonical system and consider a special ensemble formed by replicas 
of this system. The intensive properties of the system depend separately on the number of 
particles A^ on a sphere and the total area A = tcR^ of the sphere instead of being functions of 
the particle density p = N/A only. This formalism leads to the introduction of an additional 
"pressure" on top of the familiar spreading pressure p of the fluid on the curved substrate, 
which invalidates for small enough systems the Gibbs-Duhem equation p^ 128] . 

The dependence on the finite size and finite number of particles can be illustrated by 
looking at the virial coefficients of the density expansion of the spreading pressure. This 
expansion reads |1] 

/3P 



i + J2bAn,R)p'-\ (3.1) 

i=2 

where /3 = I/UbT and the first virial coefficients Bj{N, R) for a canonical system are given 



P 



by ID 



B^{N, R)={l- ^) [b^,m{R) - f 52,m(^)) ' (3-2) 

+ -^ B2,m{R) {B2,m{R) - -^B2m{R} j 



Bs{N,R)=(l-^ 



with B'-j^j{R) = dBjM{R)/dR. The Bj^m^R)'^ are the Mayer cluster integrals which are 
defined as usual, for a spherically symmetric pair interaction potential u{r) by 

B^MR) = -w^. I I dS.dS^fu (3.4) 

ZA[K) Ja{r) Ja{R) 

and 

BsMR) = -VJn^^ [ [ [ dS,dS2dSsfi2fi3f23, (3.5) 

oA[K) Ja{R) Ja(R) Ja{r) 

with fij = exp{—l3u{rij)) — 1, r^j being the geodesic distance between points i and j on 
5*2, A{R) = AttR"^, and dS is the differential area on 5*2. The first virial coefficients have 
been computed for hard disks on 5*2 (and hard spheres on S3) and the resulting approximate 
equation of state for the spreading pressure compared to simulation results. This will be 
discussed later on. Note that we have interpreted above the interaction between particles 
as acting in curved space and therefore depending on the geodesic distance on 5*2. This 
corresponds to the "curved line-of-force" case[13j. When the curved space can be embedded 
in a higher-dimensional Euclidean space (here, the sphere 5*2 in the 3-dimensional Euclidean 
space E^), one can also envisage interactions acting through the embedding Euclidean space, 
which corresponds to the "Euclidean line-of-force" case [13]. Which case arises depends on 
the physics of the system under consideration. Most existing studies use the "curved line- 
of-force" description and this is what we shall do throughout this article. 

It is worth stressing again that spherical geometry does not allow one to study finite- 
size effects at both constant curvature and constant particle density. This however may be 
viewed as a boon if one wishes to study the flat case in the thermodynamic limit. Spherical 
boundary conditions with spheres of increasing radius indeed offer an alternative to the more 
common periodic boundary conditions implemented directly in the Euclidean space. 
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C. Hyperbolic geometry 

If finiteness is the characteristic of spherical substrates, a quite different and singular 
behavior takes place in hyperbolic world. The hyperbolic plane H21 which is a homogeneous, 
simply connected two-dimensional manifold of constant negative curvature, is of infinite 
extent [291 [30] . (see Appendix A). The thermodynamic limit, at constant curvature and 
constant particle density, can therefore be taken in such a space. However, a peculiar 
feature of the hyperbolic geometry is that because of the exponential character of the metric 
at large distance, the boundary of a finite region of H2 grows as fast as the total area of this 
region when the size of the latter increases. To be more specific, consider the hyperbolic 
plane of Gaussian curvature K = —h?. The metric in polar coordinates (r, 0) is given by 

ds- = dr^ + {^^^^\c^\ (3.6) 

This form makes apparent the connection with the spherical metric of the sphere 5*2 which 
is obtained by replacing in the above expression k by iR~^ . See Appendix A. We shall refer 
to K^^ as the "radius of curvature" of the hyperbolic plane. Still in polar coordinates, the 
differential area is expressed as 

dS = ^^^^^^rfrrf0, (3.7) 

which again is the counterpart of the expression for S2 with the replacement of /t by iR~^. It 
is a simple exercise in hyperbolic geometry to show that the area of a disk of radius r is given 
by A{r) = 27rK~^(cosh(/€r) — 1) whereas its perimeter is equal to P{r) = 27r/i;~^ smh.{K,r) (see 
Appendix A). For a large radius such that nr 3> 1, one then finds as announced above that 
both A{r) and P{r) grow exponentially with the distance r, as exp(Kr). 

As a result of the above property, the boundary effects are never negligible, even in the 
thermodynamic limit! An illustration is provided by the ideal gas behavior [211 E2]- The 
canonical partition Q^{^ can be obtained from the eigenenergies of the Schrodinger equation 
for a single particle in a large domain S oi H2. This requires determining the spectrum of 
the Laplace-Beltrami operator with specified boundary conditions. For an open domain E, 
the result, obtained as an expansion in the inverse of the de Broglie thermal wavelength 



At = ^y2Trh?/{mkBT) is as follows 



,,,^, _ A(S) ^ P(S) 

XI 4Ar ' 127r 



«'^' = ir±i)i7 + 



dSK+ / dsKg 



+ 0(At), (3.8) 
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where K = —k^ is the Gaussian curvature, Kg the geodesic curvature of the boundary 
(9S, and the plus or minus signs of the second term correspond to Neumann or Dirichlet 
boundary conditions for the wavefunction on dH. The perimeter P(S) being of the same 
order as the area A(T?) for a large domain, the above expression proves the dependence of 
the ideal gas partition function on the boundary conditions in the thermodynamic limit. 

If one nonetheless insists on defining "bulk" thermodynamic properties in 7^2? a solution 
is to consider periodic boundary conditions. The procedure however is quite complex in H2 
as one must change the boundary condition when increasing the area of the primitive cell in 
which the fluid is embedded. This is detailed in Ref. [31] and is briefly explained in Appendix 
B. Imposing periodic boundary conditions amounts to replicating a primitive cell (chosen as 
a regular polygon) to realize a tiling of the whole hyperbolic plane. From a topological point 
of view, this procedure leads to a compact "quotient space" which is obtained by identifying 
in a specific manner the edges of the primitive cell by pairs. In the Euclidean case, it is well 
known that one obtains in this way for square and hexagonal periodic boundary conditions 
a one-hole torus. For the hyperbolic case, tori of genus g (i.e. with g holes) with g > 2 
are generated [Mj. Under such periodic boundary conditions, the boundary terms associated 
with 9S (and -P(S)) disappear as the manifold is now compact (boundaryless) and one is 
left with 



Qf{T) = A(S) 



1 k2 ^ 



XI 127r 



(3.9) 



which is positive and proportional to the area of the system if the de Broglie wavelength is 
small enough compared to the radius of curvature k,~^. Note that due to the Gauss-Bonnet 
theorem[7j, the area A of a primitive cell associated with a compact quotient space of genus 
g is fixed and equal to 

A = 2nK-^{g-l), (3.10) 

where, we recall, g > 2 on H2. 

This, rather long, detour via the ideal-gas limit illustrates that (i) boundary effects al- 
ways have to be considered in hyperbolic geometry, even in the thermodynamic limit, and 
(ii) "bulk" thermodynamic quantities can be defined by using periodic boundary conditions. 
Actually, a poor man's way of studying such bulk properties is to follow a procedure com- 
monly employed in analytical studies of statistical-mechanical systems defined on a so-called 
"Bethe lattice "[351 ES]. This amounts to focusing on the "deep interior" of a very large do- 
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main, far enough from the boundary and restricting in all spatial integrals that appear in 
the calculation of the fluid properties to this deep interior. In practice, one should exclude 
from computations a boundary region, which is taken as a region of width (or width whose 
ratio to the linear size of the system goes to zero when the latter goes to infinity (near the 
boundary of an otherwise very large domain), and subsequently take the thermodynamic 
limit. We shall come back to this point in the following sections. 

IV. LIQUID STATE THEORY IN SPACES OF CONSTANT NONZERO CURVA- 
TURE 

A. Statistical mechanics 

We consider a liquid at equilibrium in canonical conditions (fixed temperature T, volume 
V and number of particles A^) in a rf— dimensional Riemannian manifold. The latter is 
equipped with a metric that can be expressed in a d— dimensional set of coordinates x = 

(xi, .., Xd) as 

d 

ds^ = ^ gij(x.)dx'dx^, (4.1) 

where gtji'x.) is the metric tensor, from which one can also derive by standard differen- 
tial geometry the geodesic equation, the Levi-Civita connection and Riemann (curvature) 
tensor [71 EZl SB]- Thee infinitesimal element of volume then reads 

d 

dS = ^/\^\l[dx, (4.2) 

where ^'(x) is the determinant of the metric tensor. The geometry is assumed to be frozen, 
i.e. Qiji'x) is fixed in each point of the manifold and is not influenced by the behavior of the 
liquid. For simplicity, we restrict the presentation to a one-component atomic liquid with 
pairwise additive interactions u{r) that only depend on the geodesic distance between atoms 



computed with the metric in Eq. (4.1) ("curved line of force", see Sec. III). The canonical 
partition function is given by 

Qn{V,T) = Q^^{V,T)^^^^ (4.3) 
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where Q^^{V,T) = ^ y, , is the ideal-gas contribution aheady considered in the previous 
section and Z]^(y,T) is the configurational integral defined as 



Zn{V,T) 



V 



dSi...dS]\fexp 



N 



Yl "(^^j) 



«,i=i 



(4.4) 



with dSi given by Eq. (4.2) and the prime on the double sum indicating that the term i = j 



is excluded. As was discussed in some detail in Sec. Ill , a rigorous definition of the thermo- 
dynamic quantities may be subtle in curved space due either to the finiteness of the system 
(spherical geometry) of to the importance of the boundary conditions (hyperbolic geome- 
try). In addition to this fundamental issue, there is a practical difficulty when considering 
spaces which are not homogeneous and isotropic, for instance manifolds with spatially vary- 
ing curvature. Thermodynamics may then still be defined but the intrinsic inhomogeneity or 
anisotropy of the embedding space makes practical calculations very complex: for instance, 
the harmonic analysis (Fourier analysis, convolution theorem, Laplace eigenvalues, etc) is 
only defined on symmetric spaces, namely manifolds of constant curvature. As far as we 
know, there are virtually no existing statistical-mechanical studies of this kind. Therefore, 
in what follows, and with the goal to illustrate the infiuence of the substrate curvature on 
the behavior of the fiuid, we only consider homogeneous manifolds of constant curvature. 

From the canonical partition function, one has access to the free-energy, the energy and 
the thermodynamic pressure. The latter is defined as 

'51n(Q^\ 

/ T,Af 



f3P 



dV 



(4.5) 



As noted in Sec. |III[ the above expression corresponds to the "spreading pressure" for a 
spherical substrate. In the case of a hyperbolic geometry, we will focus on the "bulk ther- 
modynamic pressure" which is defined either through the use of periodic boundary conditions 
or by restricting the configurational integrals to the deep interior of the system, far enough 



from the boundary (see however the discussion in Sec. VB). In the canonical ensemble, one 



can also define as usual the ri— particle correlation (or distribution) functions through 



9n 11^1)1^2, •••,!■„; 



V 



N 



'N Jv 



dSn+i---dSN exp 



N 



Yl "('^^j) 



«j=i 



(4.6) 



The above expression is valid in the spherical geometry for a large enough number A^ of par- 
ticles (otherwise, there are 1/A^ corrections: see Refs[^[T3]) and in the hyperbolic geometry 
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for the "bulk" functions (see above). In both cases, the embedding space can be considered 
as homogeneous and isotropic so that gj^\ri) = 1, g]^'{ri,r2) only depends on the geodesic 
distance ri2, etc. 

Having defined the thermodynamic quantities and the correlation functions for fluids 
constrained in curved spaces (of constant curvature), one can derive the various relations 
among them. These relations take simple forms in the case where the interactions are 
pairwise additive (as considered here) and then only involve the pair correlation function. 

B. Thermodynamic pressure and equation of state 

Whereas the expression for the excess internal energy in terms of the pair correlation 
function is a straightforward extension of the standard formula in Euclidean space [SH], that 
for the thermodynamic pressure requires more caution in its derivation. We treat the case 
of positive curvature (spherical geometry) and negative curvature (hyperbolic geometry) 
separately. Indeed, in the former situation, a change of the total volume, as required from 



the thermodynamic definition in Eq. (4.5), implies a change of curvature while the total 
volume can be changed at constant curvature in the latter case. For ease of exposition, 
we deal with two-dimensional manifolds, but the reasoning is easily extended to higher 
dimensions. 

Let us consider first a spherical substrate formed by the surface 5*2 of a sphere of radius 
R. The spherical metric can be conveniently expressed in polar coordinates with two angles 
9 (colatitude) and (longitude) as 

ds^ = R^ sm{e)ddd(l), (4.7) 

or equivalently, by introducing the geodesic distance to the north pole r = R6 (0 < r < nR), 
as ds"^ = [i?sin (-^)] drdcj), which makes the comparison with the Euclidean and hyperbolic 



cases more direct. The configurational integral in Eq. (4.4) can be reexpressed as 



/"IT i*2tv /"IT i*2tt 

Zn{A,T)=R^^ / de^sm{9i) / rf0i... / d^jvsin(ejv) / rf^ivexp 
Jo Jo Jo Jo 



N 

(4.8) 






where 6ij = 6i — 6j and the total area A = AnR . By using Eq. (4.3) and Eq. (4.5) with the 
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derivative with respect to the "volume" replaced by dA = 2 (-^) dR, one finds 

(4.9) 



,P. « Az„ 



2AZjv dR 



NT 



which, after inserting Eq. (4.8 ) and the definition of the pair correlation function in Eq. (4.6 ), 
gives 

d(l)j u'iROij) % g^^\Re,, 0„ Re,, 0,), (4.10) 



where u'{r) = -^u{r). Finally, by using the homogeneity and isotropy of space, introducing 
the particle density p = N/A, and considering A^ ^ 1, one arrives at the following expression 
for the equation of state [H fT3] 



1 / dr r sin I — ] u {r)g[r), (4-11) 

2 Jo ^-R^ 



PP ^ f37ipR r\.^^„,„ fr 
P 



where g{r) = g]^ (r) is the radial distribution function depending on the geodesic distance 
r between pairs of particles. The expression can be generalized to small A^ systems and to 
high dimensions as welljU H0Hl2] . 

We now move on to the hyperbolic plane H2. The total area can be varied at constant 



(negative) curvature K = —k, but to avoid boundary problems (see Sec. Ill) in computing 



the bulk thermodynamic pressure, we use a variant of the Green-Bogoliubov method ^31 HI] . 



The derivative in Eq. (4.5) is performed via an affine transformation of the elementary area 
element, 

dS' = {l + OdS, (4.12) 



with ^ an infinitesimal parameter. Due to the form of the hyperbolic metric (see Sec. Ill and 
Appendix A), this transformation leaves the polar angle unchanged, whereas the radial 
coordinate r becomes 

^,^^^ (coshM-l) ^ (4.13) 

KSinh[Kr) 

which, contrary to the Euclidean case, amounts to a nonlinear transformation of the coor- 
dinate. We then consider the infinitesimal variation of the configurational integral which is 
generated by the above transformation. Taking advantage of the homogeneity and isotropy 
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of space (in the bulk) , one finds 

S\n{Zr,{A,T)) = NC 



(3p' 



A / dS g{r) 5u{r) 



(4.14) 



where dS = k ^ smh.{Kr)drd(j) and 6u{r) is the infinitesimal change of the pair potential 



which, by using Eq. (4.13), can be written as 

,{cosh.{Kr) 



5u{r) = ^- 



fi;sinh(Kr) 



V(r) + 0(a. 



(4.15) 



The final result reads 



P 



1 



(3np 



2k 



,.2 



dr{cosh.{Kr) — l)u'{r)g{r), 



(4.16) 



where we have taken the thermodynamic limit by letting the range of integration go to 
infinity. The above expression is the equation of state for the bulk thermodynamic pressure 
in H2. In the specific case of Coulombic systems, other types of pressures can be defined |15]: 
this will be discussed in sectioi fVBl 



Note that when the curvature goes to zero, i.e. when i? — )■ 00 in Eq. (4.11) and when 



K — )■ in Eq. (4.16), both the spherical and hyperbolic equations of state reduce to the 
Euclidean one, 

p 4 



drr'^u'{r)g{r). 



(4.17) 



However, the results on H2 and 5*2 are not simply related by replacing n hj iR ^ (as it is 



true for the metric, see Sec. Ill and Appendix A). 

To illustrate the infiuence of curvature, one may consider the virial expansion of the 



equation of state, which describes the low-density fiuid: see Eq. (3.1). The virial coefficients 
can be derived by a direct expansion of the partition function or by inserting the density 



expansion of the radial distribution function in the equation of state, Eqs. (4.11) or (4.16). 
The first correction to ideal-gas behavior is obtained by setting g{r) = e~^'"''^'^\ which leads 
after integrating by parts to 

-nR r^ 



£■ 



2I52 



dr I sin 



R 



r f r 



)(< 



-I3u{r) 



1) 



(4.18) 



and 



B. 



21^2 



- / rfrsinh(fi:r)(e-^"(")-l) 
"^ Jo 



(4.19) 
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One checks that both expressions reduce to the Euchdean formula for the second virial 
coefficient when the curvature goes to zero. By introducing as in Sec. Ill the usual Mayer 
integral with the appropriate metric, 

B2,m{R or «:-i) = -lj dSf{r) (4.20) 

with /(r) = e"''"*^'"'' — 1, it is easy to check that -B2IH2 is equal to the corresponding Mayer 



integral, whereas -B2I52 verifies Eq. (3.2) (with, here, A^ ^ 1). The difference of behavior 
between H2 and S2 reflects the property that the spherical manifold is flnite and that its 
total area can only be varied by changing the curvature, i.e. the radius of the sphere. 

C. Correlation functions and integral equations 

In liquid-state theory, it is well established that density expansions are only of limited 
use for describing the liquid phase and that a more fruitful approach consists in deriving 
approximate integral equations for the pair correlation functions [3 9J. This is based on the 
Ornstein-Zernike equation that relates the radial distribution function g{r)^ or more precisely 
the so-called total pair correlation function h{r) = g{r) — 1, to the direct correlation function 
c{r): 

h{r) = c(r) +p f dS'h{r') c(t(r, r')), (4.21) 

where S is either 5*2 and H2 (again, the formalism is easily generalized to higher-dimensional 
non-Euclidean Riemannian manifolds of constant curvature) and where t(r, r') is the modu- 
lus of the (geodesic) displacement associated with an element of the spherical or hyperbolic 



translation group. Eq. (4.21) should be considered as a "bulk" property in the case of i72, 
with all boundary effects removed when taking the thermodynamic limit. In the case of 
5*2, where flniteness is a source of additional difficulties compared to the Euclidean space, 
we consider large enough systems so that the number of particles N can be taken as a 



continuous variable and the explicit 1/A^ corrections can be neglected (see Sec. III). The 
Ornstein-Zernike equation can then be interpreted as expressing the Legendre transform 
between canonical and grand-canonical ensembles, the direct correlation function c{r) be- 
ing the second functional derivative of the grand potential with respect to local density 
fluctuations [39]. Again, this interpretation, as well as the equivalence between thermo- 
dynamic ensembles, is rigorous in the case of bulk properties in hyperbolic geometry but 
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requires some caution in spherical geometry when small systems are investigated |13t 128] . 
With this proviso in mind, one can derive the compressibility relation, 

^ = l + p f dS h{r) = \ (4.22) 

where xt is the isothermal compressibility and the second relation is a consequence of the 
Ornstein-Zernike equation. 

For homogeneous spaces of constant curvature, it is convenient to reexpress the Ornstein- 
Zernike equation by using a generalization of the Fourier transform. Harmonic analysis on 
such non-Euclidean spaces has been developed and, in a nutshell, the plane waves used in 
conventional Fourier transform for a Euclidean space are replaced by the eigenfunctions of 
the appropriate Laplace-Beltrami operator in curved space (see Appendix A). On the sphere 
5*2, one makes use of the Fourier expansion in terms of spherical harmonics expressed with 
the two angles 9 (colatitude) and (longitude). Any function f{9,(p) (recall that r = R9) 
can be written as 

oo 

f{9, ^) = J2Y1 />' O^fcK^, 0), (4.23) 

fc=0 |Z|<fc 

where 

/>, /) = r d9 smi9) [ ^ d<Pf{9, 0)y;,(^, 0) (4.24) 

Jo Jo 

and a star denotes the complex conjugate. For a function which is independent of (namely 



an isotropic function), Eq. (4.23) reduces to a Legendre transform 

oo 

f{9) = Y,f\Pk{cos{9)), (4.25) 



A;=0 

where 



/fc=f^— j / dxPk{x)f (aiccosix)) (4.26) 

and Pk{x) is the ^th Legendre polynomial. Convolutions are easily calculated so that the 
Ornstein-Zernike equation takes the form|l6] 



1 + P^^ ^hk = r^jo — . (4.27) 



It is easily checked that AtiR^Iiq = f^ dS h{r), so that the k = term of Eq. (4.26) provides 



the compressibility via Eq. (4.22) [,46]. 
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The counterpart of the Fourier transform on the hyperbohc plane H2 is the Fourier- 
Helgason transform [371 IS]- Fo^^ ^ generic function / of the polar coordinates (r, 0), it is 
defined as [57] 

f{n,t) = K-' dr c^0sinh(«:r)e^"'^P:?i/2+,,(cosh(«r))/(r,0), (4.28) 

Jo Jo 

whereas the inverse Fourier-Helgason transform reads 

/(r,0) = ^ E (-1)" / ^" tanh(7rt)e-"<^p-;)2+.t(cosh(Kr))/(n,t). (4.29) 

n=— 00 " 

For an isotropic function /(r), the Fourier-Helgason function reduces to a Mehler-Fock 
transform [48]: the dependence on n disappears and, after rewriting t = -, one has 

POO 

f{k)=2TTK~' drsmh{Kr)P_,/2+ik/.icos\i{Kr))f{r), (4.30) 

and 

/(r) = — / dkktanhl — \P^i/2+ik/.{coshiKr))fik). (4.31) 

In the above expressions, P"li/2+iti^) is a Legendre function of the first kind (conical function) 
and P-i/2+it{x) = P^i/2+iti^)- The Fourier-Helgason transform satisfies the convolution 



theorem, so that the Ornstein-Zernike equation, Eq. (4.21), can be rewritten as 



l + phk = -^^. (4.32) 

1 - PCk 



Both Eqs.(4.32) and (4.27) converge to the usual Euclidean equation with j^^^hk and h{k) 
reducing to the standard Fourier transform when the curvature goes to zero. 

One again encounters the peculiarities of curved spaces that have already been pointed 
out. We have mentioned above the difficulty associated with finiteness in the spherical 
geometry. A quite different feature occurs in hyperbolic space. One indeed finds that the 
k = component of the Fourier-Helgason/Mehler-Fock transform, 

h(k = 0) = 27TK-^ drsinh(fi:r)P_i/2(cosh(/€r))/(r), (4.33) 

Jo 

is different than the integral of h{r) over the whole space, i.e., 

h{k = 0)^2TTK-^ drsmh{Kr)f{r). (AM) 

Jo 
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Actually, h{k = 0) is always smaller than the integral. As a result, the k = component of 



Eq. (4.32) does not give the compressibility of the fluid. We shall dwell more on this issue 
when considering the liquid-gas critical point. 

The starting point of the integral equation approach to the structure and the thermody- 
namics of liquids is that the direct correlation function c{r) is simpler and shorter-ranged 
than h{r). Is is therefore a better candidate for devising approximations. Common approx- 
imations that can also be used in homogeneous curved spaces are the Percus-Yevick (PY) 
and hypernetted chain (HNC) closures, in which 

c{r)= (l+7(r))[exp(-/3M(r))-l] (PY) (4.35) 

c{r) = exp{-/3u{r) + 7(r)) - (1 + 7(r)) (HNC) (4.36) 

where 7(r) = h{r) — c{r). Once the solution of the Ornstein-Zernike equation is obtained, 
there are several routes to compute the thermodynamics from the pair correlation functions: 



the expression of the excess internal energy, the compressibility relation (Eq. (4.22)) and the 



equation of state for thermodynamic pressure (Eqs. 4.11) or (4.16)). As usual, the approxi- 
mate nature of the approach leads to thermodynamic inconstancy and the thermodynamic 
pressure obtained by thermodynamic integration of the compressibility relation 

f^ = l-- r dp'p' f dS c(r, p'), (4.37) 

P P Jo Je 

where c{r, p) is the direct correlation function at density p, differs form the "virial" expression 
in Eqs. (4.11[) or (4.16). The numerical procedures for approximate integral equations have 



been given in Ref .|16] for the sphere 5*2 and in Ref .[34j for the hyperbolic plane H2. Results 
will be discussed in Sec. |Vl 

D. Interaction potentials 

We still restrict ourselves to pairwise additive interactions and spherically symmetric po- 
tentials. If needed, these restrictions can be lifted, but it would lead here to unnecessary 
complications. We have already mentioned one subtlety concerning interactions in curved 
space: when the manifold in which the fluid is confined is embedded in higher-dimensional 
Euclidean space, distance-dependent pair potentials can either depend on the geodesic dis- 
tance ("curved line of force") or the Euclidean distance, e.g. the chord between two points 
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on a sphere ("Euclidean line of force"). To our knowledge, work on the latter type is how- 
ever very limited [TH| 1^ . Note that the alternative does not arise in hyperbolic geometry 
which is not embeddable in Euclidean space. Another subtlety in spherical geometry, and 
more generally, on compact manifold, comes with the definition of interacting pairs. Two 
particles on a sphere (or a hypersphere) are joined by two segments of a geodesic (great cir- 
cle) and can therefore interact "twice" with pair potentials associated with the two geodesic 
distances. For hard spherical particles, this clearly corresponds to the physical situation. 
For longer-range potentials and for a more formal use of "spherical boundary conditions" as 



a way to approach the thermodynamic limit in Euclidean space (see Sec. Ill and Ref. |1]) this 
"double distance convention" |H |12] is to be contrasted with the minimum image convention 
used for the "toroidal condition" corresponding to the standard implementation of periodic 
boundary conditions in Euclidean space. Note again that this feature is absent in the case 
of the hyperbolic plane H2 (although the use of periodic boundary conditions will also bring 
in compact toroidal conditions). 

Before moving on to a discussion of the nature of intermolecular forces in curved spaces, 
it is worth considering the restriction put on the form of the pair interaction potentials 
by the requirement that the various configurational integrals are well defined. In spherical 
geometry, the geodesic distance between two particles is bounded by vri? so that the interac- 
tion potentials u{r) are meaningless for vri?. (Even hard core systems make no sense if the 
hard-core diameter a in the spherical manifold is of the order of the (hyper)sphere radius 
R.) In hyperbolic geometry, a constraint arises from the exponential nature of the metric 



at long distances (see Sec. Ill and Appendix A). At large separations, the correlations die 
out and g{r) ^ 1. A typical configurational integral that appears in the calculation of the 
excess internal energy or of the pressure in H2 is then 



/•oo -I nco 

I (ir sinh(/T;r)-u(r) ~ - / dre'^^u{r), (4.38) 

Jrc ^ Jtc 



where r^. is chosen much larger than the radius of curvature k"^. The integral in Eq. (4.38) 
converges if the pair interaction decreases faster than exp(— Kr) at long distance. One there- 
fore concludes that the thermodynamic limit only exists in H2 (and higher-dimensional 
hyperbolic manifolds) if and only if the spherically symmetric pair potentials vanish ex- 
ponentially fast with distance above the radius of curvature, i.e. faster than exp(— /tr). 
For instance, this precludes using nontruncated power-law potentials such as the celebrated 
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Lennard- Jones atomic model. 



After this instructive detour, let us come back to the basics of intermolecular forces [50 
On top of the short-range repulsive interactions whose origin lies in the overlap of the outer 
electron shells of the atoms, longer-range interactions can be derived from electrostatics, 
leading to Coulombic, dipolar, etc... interactions, as well as from the multipole dispersion 
interactions between the instantaneous electric moments in one atom and those induced in 
the other [50]. Roughly speaking, all the longer-range interactions can be obtained from the 
knowledge of the Coulomb potential. In infinite Euclidean and hyperbolic space (we consider 
the two-dimensional case for simplicity), the Coulomb potential v{r) at a geodesic distance 
r of a unit point charge is obtained from the Poisson equation 

Av{r) = -2TT6^^\r,(j)), (4.39) 

where A is the Laplace- Beltrami operator in polar coordinates and S^'^\r,(f)) is the Dirac 
distribution in the appropriate manifold. The boundary condition is that the potential 



vanishes at infinity. In H2 (see Appendix A), the solution of Eq. (4.39) satisfying the 
boundary conditions is[5l] 

v{r) = - In Ttanh (^)) , (4.40) 

which behaves as 2exp(— /tr) as nr — > 00. 

In spherical geometry, the Poisson equation, has no acceptable solution as it leads to a 
singularity for particles on opposite poles. A way out is to define the Coulomb potential for 
a "pseudo-charge" corresponding to a unit positive point charge and a uniform background 
of total charge — 1 O l52] - [55] (another definition consists in taking a +1 point charge and a 
— 1 point charge located at the antipodal position [SS]). The resulting potential on a sphere 
5*2 of radius i?[5] is 

t;(r) = -ln(^sin(^)). (4.41) 

Note that the hyperbolic and spherical expressions of the Coulomb potential converge to 
the two-dimensional Euclidean logarithmic potential, — ln(r/cst), when the curvature goes 



to zero. (The extension of Eq. (4.41) to the hypersphere 5*3 is given in Ref. [51 156]).) 

The above expressions of the Coulomb potential in curved space indicate the direction that 
one should take to properly define interatomic forces in the "curved line-of-force" interpre- 
tation. In the hyperbolic geometry, for which this is crucial, one should envisage multipolar 

21 



and dispersion interactions as generated by appropriate derivatives of the Coulomb poten- 
tial. As a consequence, all such pair interactions decrease exponentially at long distance as 
exp(— Kr) or faster. Consider for instance a dipole made of two opposite charges ±g located 
at positions ri and r2. The distance separating the two charges is denoted by d. By using 
the relation between geodesic distances (see Appendix A) and taking the limit when c? — > 0, 
g — !■ oo with p = qd finite, one obtains 

pcos((h) 
smh(Kr) 

where p is the dipole moment an d(f) the angle between the dipole moment and the direction 
of the geodesic joining the position of the dipole and the center O. A crude reasoning, 
based on a one-dimensional Drude model, suggests that dispersion forces give rise to a pair 
interaction potential that goes as 

, . cosh(/tr)^ /. ,oN 

"W ■ w 4 - (4-43) 

smh(Krj^ 

Therefore, at short distance, the potential behaves as 1/r^, as expected in two-dimensional 
Euclidean space, whereas decays exponentially as exp{—2Kr) when the distance r is larger 
than the radius of curvature k~^. In Euclidean space, it is known that the simple Drude 
model captures the leading dependence on r, more sophiscated and realistic treatments only 
changing the prefactor of the interaction. 

If one is interested in studying generic properties of simple liquid models in hyperbolic 
space, with no direct connection to realistic interatomic forces, one should make sure, at 
least, that the model pair potentials are truncated beyond some cut-off so that the long- 
distance behavior remains well defined. 

E. Computer simulations 

Computer simulations play a major role in liquid-state studies in Euclidean space. Un- 
surprisingly, they have also been developed for investigating fluid and condensed phases 
in curved space. Monte Carlo simulations for short-range [13], H0VH2| [57H62] and long- 
range Coulombic-like[52], |53l ESI IM] potentials in two and three-dimensional geometry, and 
more seldom Molecular Dynamic studies in two-dimensional spherical |65l - [67] and hyperbolic 
manifolds [3TI [32| 1681470] . have been performed . Early on, numerical packing protocols have 
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also been considered in spherical [7T| [72] and hyperbolic |llj geometries. In this section, we 
do not intend to give details on the methods that have been implemented but rather to 
stress the differences with the Euclidean case and the technical difficulties encountered. 

A simplification coming with spherical backgrounds is that there is no need for periodic 
boundary conditions as the whole space is closed and finite. As already mentioned, this is 
the rationale underlying the use of "spherical boundary conditions" in order to approach 
the thermodynamic limit in Euclidean space. The price to pay is that in Monte Carlo 
algorithms one must ensure a proper sampling of the (positively) curved manifold, with e.g. 
no bias toward the poles or the equator jlO], and in Molecular Dynamic studies, one has to 
be careful with the equations of motion in curved space. The technical difficulties however 
are reasonable. 

The case of the hyperbolic plane H2 (higher-dimensional hyperbolic manifolds have not 
been studied by computer simulation), which, contrary to spherical manifolds, cannot be 
embedded in higher-dimensional Euclidean space, is much trickier. First, parallel transport 
in H2 is complex, and implementing an algorithm for solving the equations of motion either 
for hard-core particles (where collisions must be handled [STi IB2] ) or continuous interaction 
potentials (where forces must computed and vector quantities added and transported [BSl EH]) 
is by no means obvious. However, the main challenge lies in the necessary use of boundary 



conditions. As discussed in Sec. Ill and further developed in the preceding subsections IV A 



IV B IV C, one is indeed interested in the "bulk" behavior of liquids in the hyperbolic plane. 



In a simulation that is anyhow constrained to finite systems, the only way to limit the 



boundary effects, which are expected to be strong in hyperbolic geometry (see Sec. Ill), is 
to implement periodic boundary conditions pi] . This amounts to choosing a primitive cell 
(which contains the physical system) such that it can be infinitely replicated to tile the 
whole plane. To ensure smoothness and consistency, the edges of the cell must be paired 
in a specific way. In if 2, the smallest such primitive cell equipped with an edge-pairing (or 
"fundamental polygon" ) is an octagon. As the area of the primitive cell is fixed for a given 
choice of boundary condition, one must change the latter, and increase the number of edges 
of the fundamental polygon, to study finite-size effects at constant curvature and particle 
density [M]. This of course does not occur in Euclidean space where, due to the absence of 
any metric-related length scale, the system size can be changed while keeping the same type 
of boundary condition. Hyperbolic boundary conditions have been developed and discussed 
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in Ref.[3l] and implemented in Molecular Dynamics in Refs. [681470] . More details are given 
in Appendix B. 

F. Self-motion and diffusion equation 

Up to this point, we have focused on static properties. We now briefly consider aspects 
pertaining to the dynamics of fluids in curved manifolds, more specifically the self-diffusion 
of particles. In Euclidean space, the self-motion of the particles in a liquid is described 
by space-time correlation functions. The simplest one is the so-called self-intermediate 
scattering function Fs{k,t) which is the Fourier transform of the self part of the density- 
density time correlation function (or self van Hove function Gs{r,t))[39\), 

N N 

Fs{k,t) = - ^(e^'^('-^W-^W)) = - 5^(cos(H(0,t))), (4.44) 

where dj{0, t) is the distance traveled by atom j during times and t and the bracket denotes 
an equilibrium canonical average [39]. When A; — )■ and t — )■ oo, Fs{k, t) is well described by a 



2d 



(d(0,t)2) 



and in the 



Gaussian, corresponding to a truncated cumulant expansion, exp 

same long wavelength and long time limit, one expects that the self van Hove function obeys 

a diffusion equation similar to Pick's macroscopic law, which leads to Fs{k, t) = exp [—k'^Dt] 

with D the coefficient of diffusion [39] . Identification of the two expressions gives the Einstein 

relation, 

2,fl=,to{MM!>}. (4.45) 

The generalization of the self-intermediate scattering function in spherical and hyperbolic 



geometries makes use of the appropriate extensions of the Fourier transform (see Sec. IV C). 
For 5*2, one finds 

F.il^^.t) = ly IP, J CO J ^^))), (4,46) 



1 = 1 \ 



TV 

where kR G N, and on H2 

N 

N 



Fs{k,t) = ^E(^-|+^! (cosh(«:d,(0,t)))). (4.47) 

i=i 

We recall that Pknix) is a Legendre polynomial and P_i^ik{x) is a Legendre function of 

the first kind, while dj{0, t) is the geodesic distance traveled by atom j between and t and 

calculated with the appropriate metric. 
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One may then wonder what is the small k and long time limit of the intermediate scat- 
tering function in non-Euclidean space. One expects that it converges to a result predicted 
by the associated diffusion equation, 



dt 



DAGs{r,t), (4.48) 



with A the Laplace-Beltrami operator acting in the curved manifold. (Actually, by isotropy 
one expects Gs to only depend on r so that only the "radial" part of the operator is required.) 
We keep focusing on the two-dimensional manifolds 5*2 and H2 for simplicity. The diffusion 
problem has been solved in both cases with the solution in 5*2 given by [73] 

Gf='(r.«) = E (^) sm (1) ft (cos Q) ,^^^, (4,49) 

fc=0 ^ ^ 

and that in if 2 by [71] 

-iiiD« /.oo -J^ 

f 4: I IIP ADt 

G[^^\r,t) = ^^ ^sinh(fi;r) / dy ^ (4.50) 

2v27r(Dt)2 Jr A/cosh(Ky) — cosh(fi;r) 



One can check that when the curvature goes to zero, i.e. R — )■ +00 in Eq. (4.49) and k — )■ 



in Eq. (4.50), both expressions reduce to the standard two-dimensional Euclidean formula. 



/n \ e 4-Dt 

«■ "('■•«)= 4^- (^'^^ 

Generalizations to the three-dimensional manifolds 5*3 and H^ are given in Ref. [731 EH] and 
Ref. [7^, respectively. 

Diffusion on a sphere (or a hypersphere) is such that the particles can never escape 
from the closed space they are embedded in and, as a result, the mean square displacement 
((i(0,t)^) goes to a constant value when t — t- +00. Procedures to fit simulation data for 
extracting a self diffusion coefficient from the ((i(0,t)^) have been given in Refs. [73] [75] . Dif- 
fusion on Riemann manifolds (with a non constant curvature) has been recently investigated 
in the short-time limit [76]. For the hyperbolic plane, the diffusion motion at long times is 
characterized by two regimes: when the mean traveled distance is large, but smaller than 
the radius of curvature k~^, one has an ordinary diffusion regime with ((i(0,t)^) ~ Dt and 
{d{0,t)) ~ vDt; on the other hand, when the mean traveled distance is large compared to 
K-\ one finds a "ballistic" diffusion regime with {d{0,tf) ~ {Dtf and (ci(0,t)) ~ DtfMiW7\. 
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FIG. 1: Log-log plot of the mean absolute displacement {d{0,t)) in units of k ^ and of 



{f{d{0,t))) (see Eq. (4.52)). The dashed line has a slope equal to 1 and the dotted-dashed, 



a slope of 1/2. To good approximation, {f{d{0,t))) has a linear time dependence at all 

times, which corresponds to a diffusive regime in the hyperbolic sense. The parameters are 

pa"^ = 0.91, T = 2.17, and na = 0.2. The system with octagonal periodic boundary 

conditions comprises 287 atoms. 



It has been shown in Ref . [77] that the self-diffusion coefficient can be conveniently extracted 
from simulation data of the function -^ '^j=i{f{dj{0,t))) where 



fir) 



ft-" In 



1 + cosh.{K,r] 



(4.52) 



Indeed, when calculated with the solution of the hyperbolic diffusion equation, this function 
behaves as Dt for both the ordinary and ballistic diffusion regimes. In Fig. [T| we show the 
Molecular Dynamics simulation result for a fluid with truncated Lennard- Jones interactions 
on H2 (see below) [77j at a rather high temperature. It is clearly found that the atomic motion 
becomes diffusive at long enough time, and the self-diffusion coefficient can be determined 
without having to worry about the crossover between ordinary and ballistic diffusion. 
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V. CURVATURE EFFECT ON THE THERMODYNAMICS AND STRUCTURE 
OF SIMPLE FLUIDS 



A. Short-range interaction potentials 

In order to illustrate the influence of a nonzero curvature on the properties of a fluid 
embedded in a non-Euclidean space, we begin with the two dimensional hard-disk fluid on a 
sphere 5*2 of radius R and a hyperbolic plane H2 of radius of curvature n~^. (We recall that 
the Gaussian curvature K is constant in both cases and equal to 1/ R^ for 6*2 and — k~^ for 
H2-) The "virial" equation of state giving the thermodynamic pressure, i.e. the spreading 
pressure for S2 and the bulk pressure for H2 (see preceding sections), is obtained from Eqs. 



(4.11) and (4.16) after using the fact that for hard disks of diameter a in curved space, 



dr 



exp(-/3M(r)) = 5{r - a). For ^aH US 



and for H2 



P 
P 



1 + p- 



7ri?2 



52 



Hi 



£1 ^^" ill) ^^^ 



n 



1 + Ptt^ (cosh {k(t) - 1) g{a^ 



(5.1) 



(5.2) 



where in both cases g^c^) is the radial distribution function at contact. 



In the low density limit, g{cr^) goes to 1 and the second viral coefficient (see also Eqs. (4.18) 



and (4.19)) reads 



B^ 



2I52 



2 [r)'''{^ 



and 



B2\ 



Hi 



-K 



K^ 



(cosh (ko") 



(5.3) 



(5.4) 



where we have considered A^ ^ 1 for 5*2 and Eq. (5.3) is only defined for a < ttR. For H2 



even if nothing a priori prevents one from taking the particle diameter much larger than n ^, 



the discussion in Sec. IV D indicates that the physical interactions decrease exponentially 
with distance beyond k"^. Particle diameters a ^ k^^, which lead to an exponentially 
growing -B2IS , have therefore no physical significance. Higher-order virial coefficients have 
been calculated for 5*2 [U [13] and for if 2 [31]. (The three-dimensional spherical case S3 has 
been considered in Refs. |iOl HH 159].) 

As already noted, there is no simple symmetry k, 4-7- iR^^ between the H2 and 5*2 pressure 
expressions, due to the finite size of the sphere whose total area can only be changed by 
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varying the curvature. For instance, when the curvature goes to zero in 5*2 and H2., the 
equation of state and virial coefficients approach the Euchdean result, albeit in a non- 
symmetric way. As an illustration. 



B^ 



2|S2 



B^ 



2Ie2 



1 - -Ka'' 
6 



O(irV^) 



B. 



Bo 



2 4^ 



1 + —Ka' + 0{K'a 



(5.5) 



(5.6) 



where B2\^ = 7ia'^/2. The asymmetry between 5*2 and H2 is not compensated if one express 
the virial expression in terms of the packing fraction rj (or rather "surface coverage " in two 
dimensions) instead of the density p. Indeed, 



so that to ffist order, 



P 



v\h,= 
= 1 + v 



27tR 

2n 



'''{^))P^ 



K 



2vC«sh(f) 



1 P, 



and 



52 

pp_ 

p 



/ M sin (^) COS (^) 
\2r) 1-cosf^) ^ ^^ 



sinh(^) ^^ 2^ 

cosh(f)-l ^'^ 



(5.7) 
(5.8) 

(5.9) 
(5.10) 



To go beyond the low-density expansion, several routes have been taken for exploring the 
hard-disk fluid in spherical and hyperbolic geometries. One can build approximate equations 
of state from scaled particle theory [IHl [TS], free area theory [211 E2], or by constructing 
rational approximants from the ffist virial coefficients [31], with possible additional input 
from some close packing density [131 EHl CS]- One can also compute the radial distribution 
function g{r), either from approximate self-consistent integral equations[lH HH] of from 
simulations [m [311 [321 [591 [HO]- The thermodynamic pressure can be obtained through the 
"virial", Eqs.( 5.l| ) or (5.2), or the compressibility, Eq. (4.37), relations. 

As in Euclidean space, all approximate methods, including integral equations, are lim- 
ited to moderate densities, i.e. densities that are less than the ordering transition in the 
Euclidean plane. The exact location, the nature and the order of the latter are known to be 



strongly system-size dependent (see Sec. VI), so we typically mean here densities pcr^ < 0.8 
and packing fractions r] < 0.6. In Fig. [2| we compare for the case of the hyperbolic plane 
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FIG. 2: Equation of state of the hard-disk fluid in H^ versus packing fraction 77 for a quite 

strong curvature kg ^ 1.06. Comparison between the PY results (compressibility route 

(dotted line) and virial route (full line)) for kg = 1.06 and the simulation results of 

Ref. [31] for na = 1.060 (green squares) and na = 1.062 (blue diamonds). The 

dashed-dotted line is the prediction of Ref. [79] and the dashed line corresponds to the 

equation of state truncated at the second virial coefficient. The disagreement between 

predictions and simulation data is significant alrealdy for rj ^ 0.2. 



with a quite strong curvature na ~ 1.06 the equations of state obtained by simulation[ 
and predicted by several approximate methods (approximate equations of state [211 E21 ES] 
and results from the Percus-Yevick integral equation with both the virial and the compress- 
ibility routes [H] As in the Euclidean case, the compressibility route gives a higher bulk 
thermodynamic pressure than the virial route. An empirical recipe that has proven very 
accurate in Euclidean space ^39j is to consider a linear combination 

2Pe + Pv 



P 



(5.11) 



where Pc and Py are the "compressibility" and "virial" pressures. In Fig. |2| one can see 
that all predictions, except for the free-area method (which for this reason not reproduced 
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FIG. 3: Equation of state of the hard-disk fluid in 5*2 and in H2 for the same magnitude of 



the curvature hia 



R 



0.45. The symbols denote simulation results in 6*2 (squares) and 



in E2 (circles). The full curve corresponds to the PY approximation in H2- The PY result 
in 5*2 [IS] is virtually indistinguishable from the simulation results (on this scale). 



here), approach at low density the exact expression truncated at the second virial coefficient, 



Eq. (5.9), whereas their validity rapidly deteriorates for even moderate coverage. One should 

have been carried out with 



however keep in mind that the corresponding simulations [ 
a very small number of atoms, always less than 10, and that the curvature is quite strong. 
For the spherical substrate 5*2, most of the existing data concern the approach to the 
Euclidean plane by following the equation of state (3P/p versus 77 or pa"^ at constant number 
A^ of atoms (which implies that the curvature varies with density, as pa^ = (£:) (^) ). We 
nonetheless display in Fig. |3] a comparison between the equations of state obtained from the 
the Percus-Yevick integral equation via the compressibility route for the hard-disk fluid in 
5*9 |46j and if2|ll] for the same "radius of curvature" a/R = na = 0.45, as well as simulation 
results in S2 and E2. Despite possible limitations of the Percus-Yevick approximation, it 
appears that the results in 5*2 almost coincide with those in the Euclidean E2 whereas those 
in H2 are significantly different for pa^ > 0.55. 

30 



2.5 



S: 2 

CO. 



1.5 



0.2 



0.4 



0.6 0.8 

Ka,c/R 



1.2 



FIG. 4: Compressibility factor /3P/p of the hard-disk fluid in 5*2 and H2 as a function of 

the curvature parameter, a/R or kct. For 5*2, are plotted the result from the PY 

approximation for pa'^ = 0.8 116] (full line, where for convenience [PP/p — 5) is plotted) and 

the simulation data for 77 = 0.4 (which corresponds to pa"^ varying between 0.385 and 0.392 

J13j (symbols). For H2, the two curves are obtained from the PY approximation with 

pa"^ = 0.2 (dotted line) and pa"^ = 0.4 (dashed line). 



In addition, to illustrate the influence of curvature, i.e. both its magnitude and its sign, 
on the thermodynamic properties of a fluid, we have collected in Fig. |4] simulation data and 
integral-equation results for the compressibility factor pP/p of hard disks on 5*2 and H2, 
which we plot as a function of increasing curvature, \K\~^^'^cr, for different values of pa^. 
Whereas the pressure increases with curvature in 5*2 for a given density pa"^ (or packing 
fraction rj), the opposite trend is observed in H2. This effect can already be seen as very 
small densities when considering only the correction to the ideal-gas behavior due to the 
second virial coefficient (see above). 

The main advantage of computing the radial distribution function g{r) is that it also 
provides some information on local structural order in the fluid. We focus here on the fluid 
or liquid regime at densities and temperatures (for continuous interaction potentials) such 
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that the system in flat space has not undergone its ordering transition (and if attractive 
interactions are present, away from the gas-hquid critical point). Higher densities and lower 



temperatures, which correspond to interesting new physics, will be considered in Sees VI and 



VII, and critical behavior in Sec. VC In the fluid/liquid regime at "moderate" densities and 
temperatures (the qualifler "moderate" having the above discussed meaning), it is found that 
curvature has only a weak influence, at least when the curvature is not too large, i.e. when 



-R'Icr is signiflcantly less than 1. This can already be inferred from the thermodynamic 
data (see above). We illustrate the effect of a negative curvature on the g{r) of the hard- 
disk fluid, as obtained from the Percus-Yevick equation, in Fig. [5144J. For a packing fraction 
rj = 0.55, it is found that the curves are essentially superimposable on the Euclidean one for 
a range of curvature parameter na between and 0.5. Only when going to higher curvatures, 
e.g. KO" = 1.5, can one distinguish the influence of the curvature in the structure. (Note 
that on the sphere, strong curvatures with o" > i? correspond to small systems for which 
flnite-size and discreteness effects are very important|13].) 

The same conclusion about the curvature effect has been reached for a truncated Lennard- 
Jones model in H2, with no detectable changes from flat-space results up to na ~ 0.5| 
The model is deflned by a pair potential 



u{r) = 4e 



cr\l2 /cr 






Uc (5.12) 



that is truncated for r > Vc with Vc = 2.5cr; Uc is the shift obtained from the relation 
u{rc) = 0. In this study, a large range of temperature, density and curvature parameter has 
been investigated through the solution of the Percus-Yevick and HNC integral equations. 
As far as we know, this is also the only work in which the validity of such integral equations 
in curved space has been tested. As shown in Fig. 7, of Ref . [H] , the predictions compare 
well with Molecular Dynamics simulation data with a slight advantage to the Percus-Yevick 
approximation. 

B. Coulombic systems 



We have briefly reviewed in Sec. |IVD| the way to deflne the Coulomb potential in curved 
space (following a "curved line of force" interpretation) via the solution of the relevant Pois- 
son equation. We have also stressed the difference between hyperbolic geometry where space 
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FIG. 5: Radial distribution function g{r) of the hard-disk fluid in H2, as obtained from the 

PY equation, at a packing fraction rj = 0.55 for various values of the curvature: 

Ka = 0, 0.15, 0.5, 1.5. Only for k = 1.5, namely a radius of curvature smaller than the 

particle diameter, does the structure display significant deviations from the Euclidean case. 



can be infinite and spherical geometry in which the finiteness of space requires the introduc- 
tion of pseudo-charges in order to properly define the Coulomb potential. The main advan- 
tage of spherical geometry being the absence of boundary, "spherical" and "hyperspherical" 
boundary conditions have been widely used by Caillol, Levesque and coworkers [5|, El E2}- 
\55\ E21 El] to approach the thermodynamic limit of Coulombic systems in Euclidean space 
(as the radius of the sphere or the hypersphere goes to infinity). Multipolar interactions 
deriving from the Coulomb potential have also been considered [53]. Here, as in most of the 
article, we focus on two-dimensional manifolds of constant positive {S2) or negative (-^2) 
curvature and we consider systems of charges interacting through the Coulomb potential. 



'^i^)\s2 = ~'?i^2 In f sin I — 1 



\2RJ 



uii")\s2 = -?i?2ln (tanh (yj j 



(5.13) 
(5.14) 
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for two charges qi,q2 (or pseudo-charges in S2) separated by a geodesic distance r (with 
r < 71 R in 6*2). 

The two simplest Coulombic systems are the one-component plasma (OCP) and the 
two-component plasma (TCP). The former is a monodisperse system of point charges of 
equal sign and magnitude q embedded in a charged uniform background that maintains 
global electroneutrality. Conventionally, one uses the coupling constant F = /3g^ as control 
parameter. At small and intermediate coupling, F = 0(1) or less, the system is in a fluid 
phase whereas at zero temperature (F — ;■ +00), its ground state in the Euclidean plane is a 
hexagonal crystal forming a triangular lattice. The questions which remain debated [67j are 
(i) whether there is a freezing transition at some large but finite coupling and (ii) if indeed 
freezing occurs, whether the transition is a first-order one or proceeds via two continuous 
transitions separated by a hexatic phase as predicted by the KTNHY theory |9l |80l [81]. As 



will be more extensively discussed in Sees VI and VII , placing the system in curved geometry 
"frustrates" hexagonal order by forcing in an irreductible number of topological defects. We 
just mention here that it has then been argued that approaching the thermodynamic limit 
of the OCP in the Euclidean plane via spherical boundary conditions avoids the artifacts 
generated by periodic boundary conditions on the ordering behavior [67j. 

The TCP consists of a binary mixture of oppositely charged point particles (with charges 
±g) in equal concentration. No neutralizating background is needed. In the Euclidean 
plane, the system is stable in a conducting phase up to a coupling F = 2[H21 ES]. Above this 
value, a collapse of pairs of opposite charges occurs. This can be regularized by introducing 
an additional hard-core repulsion (the system is then also known as the restricted primitive 
model). The system then undergoes a continuous Kosterlitz-Thouless transition at F = 4 to a 
dielectric phase in which all charges are bounded into dipoles in the limit of an infinitesimal 
particle density. For larger densities, the critical temperature of the Kosterlitz-Thouless 
transition decreases and terminates in a first-order transition slightly above the gas-liquid 
critical point. The location of the latter remains quite difficult to obtain precisely, as both 
the temperature and the density are very small [84j- 

In the following, we restrict ourselves to the OCP and TCP in their fluid, conducting 
phase at small and intermediate coupling in spherical (6*2) and hyperbolic geometry (-^2)- 
Several properties found in the Euclidean plane carry over to 5*2 and H2: (i) the small 
coupling/high temperature regime is (asymptotically) described by the linearized Debye- 
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Hiickel approximation, (ii) exact solutions for the thermodynamics and the structure are 
obtained for the special value of the coupling constant F = 2, and (iii) in the conducting 
phase, exact sum rules generalizing the Stillinger-Lovett relations are satisfied. 

We begin with the OCP. Assuming a perfect compensation between the point charges 
and the background amounts to replacing g{r) by the pair correlation function of the point 



charges, h{r) = g{r) — 1, in the equation of state for the thermodynamic pressure, Eqs.(4.11 ) 



and (4.16). After inserting Eqs.(5.13) and (5.14), one finds 

S2 



P 
P 



1+P 



n 



r 



4 



drr (1 + cos 



R 



1+P 



H2 



ttT 



dr 



cosh(Kr) 



which both converge to the Euclidean formula 

/3P 
P 



1 + P 



ttE 



sinh(/T;r) 



drr h{r) 



h{r) 



h{r) 



(5.15) 



(5.16) 



(5.17) 



E2 \ ^ / JO 

when the curvature goes to zero. It should be stressed that defining a pressure for the 
OCP is far from trivial as both the way one treats the uniform background and the choice 
of embedding conditions may matter. This has been carefully discussed in the case of the 
hyperbolic plane H2 for which no less than 5 different pressures have been considered |45]. 
The above defined "bulk thermodynamic pressure" reduces to the standard thermodynamic 
pressure (then also equal to the "Maxwell pressure") with vanishing curvature. 

It has been shown for S'2[85j and iif2[5l] that the OCP in the conducting phase satisfies 
generalized Stillinger-Lovett sum rules [86]. 



pi dSh{r) = -l, 

'si or H2 



which expresses the strictly enforced electroneutrality, and 



ttR'^P^T f dSh{r) ( 

JS2 ^ 



1 — COS 



R 



47rp2r 



H2 



dS h{r) In ( cosh ( — j 



(5.18) 



(5.19) 



(5.20) 



which both express the screening property, with ^j2'KpT the inverse of the Debye screening 
length. Note that the sum rules in the case of H2 should be interpreted as valid in the bulk 
(see Sees, [m] and llV] and Ref . [87] ) . 
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In the high-temperature or small-couphng hmit, the hnearized Debye-Hiickel approxi- 
mation becomes asymptotically exact and analytical expressions for the pair correlation 
function can be derived. The direct correlation function c{r) becomes equal to —f3u{r), so 



that by using the Ornstein-Zernike relation, Eq. 4.32, one obtains in 5*2 



hiO) 



2{x' + 



-. oo 

' fc=l 



2k + I 



X^ + k{k + l) 



PkicosiO)) . 



(5.21) 



(No explicit expression of h{6) has been obtained, contrary to the case of the hypersphere[5].) 
lnH2, one finds [EH EZl 

h{r) = -rg^(cosh (/tr)) (5.22) 



where 



1 



u 



+ 



27TpT 



(5.23) 



2 V 4 K- 

As anticipated, one can show that the above Debye-Hiickel expressions satisfy the gen- 



eralized Stillinger-Lovett sum rules |51j. Eqs.(5.21) and (5.22) can then be inserted in Eqs 



(5.15) and (5.16) to derive the equation of statejH]. For illustration, we consider the F — )■ 



limit in the hyperbolic plane. At fixed curvature, Eq. (5.23) shows that z/ — )■ and, by using 

2 

^, one finally obtains 



Qo(a:) 



the definition Qo{x) = | In (jzf ) and the property that f^ dx ^ 
that 



P 



Hi 



3f2 



pn-^T 
12^ 



(5.24) 



as F — 7- 0. On the other hand, at fixed (small) F and k — )■ 0, one finds that the Debye-Hiickel 
expression in H2 leads to 

1 - -, (5.25) 



P 



H2 



F 

4' 



which is the exact result for the thermodynamic pressure. It is also worth noting that in 
H21 the exponential decay of the Coulomb potential restores an analytical virial expansion 
in the Debye-Hiickel limit for a finite radius of curvature k^^ whereas in Euclidean spaces, 
the pressure is independent of density in E2 and has a nonanalytical behavior in E^ [88] . 

As mentioned above, for the specific value of the coupling constant F = 2, exact analytical 
expressions of the thermodynamic quantities and the pair correlation function h{r) can be 
obtained in Euclidean [89], spherical [851 [90], ^'^^ hyperbolic geometries [5T| l9T] (as well as on 
a "Flamm's paraboloid "|92j). For spherical and hyperbolic geometries, the pair correlation 
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function is given by 



and 



h{r)\ 



■l + cos(i) 



iTvpR^ 



S2 



-(cosh(|:) 



Sivp 



(5.26) 



(5.27) 



From the above expressions, one can calculate the thermodynamic pressure with the follow- 
ing results: 

RP 1 

= o (5-28) 



P 



52 



and 



P 



271 p + K^ 

H2 47rp + k2 



(5.29) 



Note that the low-density limit of Eq. (5.29) coincides, as it should, with Eq. (5.24) when 
r = 2. When the curvature goes to zero, Eq. (5.29) reduces to the Euclidean result, — = 



|, whereas the spherical expression always coincides with it, regardless of the (positive) 
curvature. An interesting observation concerning the hyperbolic plane is that if one takes 
the limit of an infinite (negative) curvature, k, — )■ oo, the pressure goes to the ideal-gas limit. 



This result is expected to be quite general (see Eq. (5.28) and Sec. VB): the influence of the 



interactions appears to vanish in the limit of large negative curvature, if however this limit 
makes any sense. 

Having discussed the OCP in two-dimensional curved manifolds (more details can be 
found in Refs.[Ml ESI ED] for ^2 and Refs.|13 iSl ED EZl E] for H2), we now move on to 
the case of the TCP model. In this system, there is no background and one has to consider 

the correlation functions between pairs of equal charges, /i++(r) = h (r), and of opposite 

charges, /i+_(r) = h |_(r). Electroneutrality implies 

p j dS {h^+{r) - h+4r)) = -2, (5.30) 

where p = 2p+ = 2p_ is the total particle density. Higher-order Stillinger-Lovett types of 
sum rules are also satisfied by /i++(r) and /i+_(r) and the thermodynamic pressure from 
Eqs.( 4.1l| ) and ( |4.16[ ) is now given by equations similar to Eqs. (5.15) and (5.16) with p 
replaced by p/2 and h{r) by (/i++(r) — /i+_(r)). 

In the small coupling limit, the Debye-Hiickel approximation gives the first nontrivial 
term in F. As the linearized Poisson equation is similar to that of the OCP, the pair 
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correlation functions are given by /i++(r) = — /i_|__(r) with /i_|__|_(r) + hj^^{r) equal to the 



solution in Eqs.(5.26) and (5.27). Consequently, the equation of state as the same of the 
OCP given above. Exact analytical results are also obtained for the special coupling value 
of r = 2. However, the calculations become quite involved and tedious. Details concerning 
the hyperbolic case can be found Refs.fSTl 1871 191] . 

C. Liquid-gas critical behavior 

Critical behavior implies fluctuations on all spatial scales and a diverging correlation 
length. Strictly speaking, critical phenomena cannot exist in finite systems and therefore 
not for systems confined to spherical substrates. (As stressed before, spherical geometry 
can however be used as a trick to approach the thermodynamic limit in Euclidean space 
and to study finite-size scaling by decreasing the curvature and concomitantly increasing 
the system size [62].) On the other hand, bona fide critical points can be present in systems 
embedded in hyperbolic geometry as space can then be infinite. We consider in this section 
the gas-liquid critical behavior of atomic fluids in the hyperbolic plane H2. 

Drastic changes in critical behavior are expected in hyperbolic geometry. This is known 
from field theoretic studies: a negative curvature acts as an "infrared regulator" , suppressing 
fluctuations on wavelengths larger than the radius of curvature |93]. As a consequence, 
critical behavior in statistical systems on a hyperbolic manifold is expected to be mean-field- 
like, with classical values of the critical exponents (r/ = 0, 7 = 1, etc) pH - [^ . This mean-field 



character can be understood by going back to the analysis of section IV C concerning the 



correlation functions and the compressibility. The latter is given by an integral over space. 



Eq. (4.22), which in the case of H2 is expressed as 



00 



— — = 1 + p — / dr sinh(fi;r)/i(r). (5.31) 

P ^ Jo 

For long distances, sinh(/€r) ~ (1/2) exp(/€r), and, as already pointed out in the section 

concerning pair interactions, the integral over the whole space H2 is bounded if and only 

if h{r) decreases at large r faster than exp(— Kr). Assume for instance that h{r) decreases 



exponentially as exp(— r/,^) with ^ the correlation length. One then finds that Eq. (5.31) 
becomes 
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which diverges for k^ — )■ 1 (and is infinite for k^ > 1). The compressibihty can therefore 
diverge with a finite correlation length ^ = k^^ (this is the meaning of the "infrared regu- 
lator" discussed above). The critical point being then characterized by a finite correlation 
length, one expects ^ to be a regular function of the control parameters T and p even in the 
vicinity of {Tc,pc), 

e(r, p) = K-' - A{T - TJ - B{p - p,f + ... (5.33) 

with A, i? > 0. As a result, along the critical isochore, 

XT~|T-T,|-\ T^T+ (5.34) 

and along the critical isotherm 

XT-- {p- PcY^, P^ Pc (5.35) 

which gives the classical values of critical exponents, 7 = 1, 5 = 3, and rj = 0. Note also 
that, at the critical point, the Fourier-Helgason transform of the pair correlation function 



h{k) (see Eq. (4.30)), is finite for A; = 0, which is a dramatic illustration of Eq. (4.34), and 
has a regular expansion in {k/n) when (fc/zt) — > 0. 

We stress that the above reasoning concerns the "bulk" behavior in H2, that which is ob- 
tained by taking the thermodynamic limit either with a succession of periodic boundary con- 



ditions or by removing the boundary regions (see Sec. Ill and IV). The exponentially grow- 
ing character of the hyperbolic metric otherwise induces possible "boundary transitions" [95] 
which we do not consider here. 

The gas-liquid critical point of the truncated Lennard- Jones fiuid in H2 has been studied 



through the integral equation approach[44j that has been described in Sec. IV C The results 
are illustrated in Figs|6]and[7]for the Percus-Yevick closure and a curvature parameter na = 
0.5. Fig. [6] displays the inverse compressibility along different isotherms as one approaches 
the critical one. It is well known that approximate integral equations such as the Percus- 
Yevick one fail to capture nonclassical critical behavior when it is present, as in the Euclidean 
plane. The relevant conclusions that can be drawn from the Percus-Yevick equation study is 
therefore not about critical exponents, but rather concerning the above discussed scenario of 
an exponentially decreasing pair correlation function with a correlation length approaching 
the radius of curvature k,~^ from below. A clear confirmation of the scenario can be seen in 



Fig. [7|where both h{r) and exp(Kr)/i(r) are plotted in the vicinity of the critical point. Note 
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FIG. 6: Inverse of the isothermal compressibility versus packing fraction t] for the 
truncated Lennard- Jones fluid in H2 within the PY approximation near but above the 

critical temperature: from top to bottom , 
T = 0.640, 0.610, 0.580, 0.560, 0.530, 0.515, 0.490, 0.472 (in the usual reduced units). The 

curvature parameter is na = 0.5. 



that the gas-liquid critical point for na = 0.5 is found at a lower temperature than that in 
the Euclidean plane. 

An interesting question concerning the critical point in the presence of negative curvature 
is whether the location of the point in the (T, p) diagram goes continuously to the value in 
the Euclidean plane or goes to a lower, possibly zero-temperature point as the curvature is 
reduced to zero. The former possibility is natural if one thinks of the curvature as introduc- 
ing a mere finite-size cut-off. However, the latter has been conjectured by Angles d'Auriac 
et a/j95] in their study of the Ising model on hyperbolic lattices; these authors have also 
suggested the existence of a crossover line, emanating from the critical point in Euclidean 
space and extending to small nonzero curvature below which the pair correlation function 
has an algebraic decay for r < k,~^ and an exponential decay at long distance. Approxi- 
mate integral equations presumably predict a continuous behavior, and no crossover is seen 
around the location of the flat space critical point, but this may come from their intrinsic 
limitation. We have undertaken preliminary Molecular Dynamics simulations of the trun- 
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FIG. 7: (a) Log-linear plot of the pair correlation function h{r) = g{r) — 1 of the truncated 

Lennard- Jones fluid in if2, as obtained in the PY approximation for T = 0.472. The 
behavior as a function of r] is nonmonotonic and the slowest decay is for r] = 0.118, which 

corresponds to the maximum of the compressibility. The dotted line is exp(— Kr). (b) 
Same data multiplied by exp(/€r): a convergence towards a constant plateau at large r is 
clearly visible as one approaches the critical point. The curvature parameter is na = 0.5. 



cated Lennard- Jones model in H2 to try to locate its gas-liquid critical point and possible 
crossover behavior. Fig. Is] displays a three-dimensional plot of the integral over space of h{r) 
in the (T, 0) plane for a curvature parameter na = 0.2. Periodic boundary conditions with a 
fundamental polygon of 14 edges (see Appendix B and Ref . j34] ) are used and the system size 
is unfortunately rather small (from 8 to 320 atoms). One observes that the integral indeed 
starts to rise very steeply as one reaches the vicinity of the Euclidean-space critical point 
and saturates at lower T and p due to system-size limitations. However, without studying 
finite-size effects by changing the fundamental polygon of the periodic boundary conditions, 
no clear-cut conclusion can be reached. Further work is needed to elucidate this question. 
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FIG. 8: MD simulation of the critical behavior of the truncated Lennard- Jones model in 

H2 for KO" = 0.2. Three-dimensional plot of the integral of h{r) over the system size as a 

function of T and 77. The critical point of the model in E2 obtained by Monte Carlo 

simulation [97j is around Tc = 0.472 and r] = 0.33. 



VI. FREEZING, JAMMING AND THE GLASS TRANSITION 

A. Geometric frustration 

Starting with the work of Frank[8], BernalpS] and others in the fifties, a whole line of 
research has developed, trying to understand liquids, glasses and amorphous packings, i.e. 
systems with no apparent structural long-range order, from a geometric point of view. An 
explanation for the avoidance of crystal formation when cooling a liquid, with the resulting 
glass formation, as well as for the structure of glasses and amorphous materials has been put 
forward in terms of "geometric frustration" [7] . The latter describes a competition between 
a short-range tendency for the extension of a locally preferred order and global constraints 
that preclude tiling of the whole space by a periodic repetition of the local structure. A 
prototypical and well documented example is that of icosahedral order in three-dimensional 
systems in which particles interact through spherically symmetric pair potentials: despite 
being locally more favorable, icosahedral order built from tetrahedral units cannot freely 
propagate in space to give rise to long-range crystalline order. Icosahedral order therefore 
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must come with topological defects. A step forward in the geometric description has been 
provided by the "curved space approach" developed by Kleman, Sadoc, Mosseri on the one 
hand[71 dU] and by Nelson, Sethna, and coworkers on the other P EH]. In the late seventies, 
Kleman and Sadoc [10] realized that perfect icosahedra could tile three-dimensional space 
if the metric of the latter were modified to introduce a constant positive curvature. On 
the hypersphere S-^ with a radius R equal to the golden number ( ^^^ ) times the particle 
size, 120 particles form a perfect icosahedral tiling known as the {3, 3, 5} polytope[7j. Such 
an unfrustrated crystallization is then expected to take place more easily, i.e. at higher 
temperature, than crystallization in Euclidean space which, due to frustration of icosahedral 
order, must involve a reorganization of the local order and a different, hexagonal close- 
packed or face-centered cubic, long-range order |571 EH]. The template or "ideal order" can 
then be used to describe real physical systems in Euclidean space, mostly metallic glasses, 
provided that one can described how it evolves when "flattening" space: roughly speaking, 
one observes in Euclidean space the remains of ideal icosahedral order that is broken up by 
the necessary appearance of topological defects, essentially disclination lines, which in some 
sense carry the curvature needed for compensating the flattening of the template. The same 
approach can be taken for a variety of local orders that are frustrated in Euclidean space 
but can tile space in curved spherical or hyperbolic geometries [7] . A theory of the glass 
transition has also been developed along these lines |100] llOlj . 

Quite generally, curvature can be used as a tool to either frustrate or generate long-range 
order. In the above discussed example of icosahedral order, curvature allows one to build 
a perfect tetrahedral/icosahedral tiling in the form of a polytope in S^. Generically, such 
polytopes can be found for special values of the radius of curvature expressed in units of 
the polytope edge: there are a finite number of them in spherical geometries and an infinite 
number in hyperbolic geometries (one could also consider geometries with spatially varying 
curvature) 

Consider now two-dimensional space and fluids of particles interacting via spherically 
symmetric pair potentials. Such systems of disks on the Euclidean plane are not subject 
to geometric frustration: the locally preferred structure is a regular hexagon, with one 
atom at the center and 6 neighbors at the vertices, and this structure can be periodically 
repeated in space to form a triangular lattice. The system crystallizes extremely easily, 
either through a sequence of two continuous transitions with an intermediate hexatic phase 
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as in the KTNHY[9l EOj |8T] or through a weak first-order transition, and glass formation 
never occurs. Curving space then frustrates hexagonal order and forces in topological defects 
which are point-like in two dimensions. (The topic of hexagonal order and associated defects 



will be considered in more detail in Sec. VII ) Again, for specific values of the radius of 



curvature compared to the particle size, the spherical and hyperbolic manifolds (5*2 and 
if 2, respectively) allow crystalline-like tilings which are conventionally represented by the 
Schlafii notation {p, q} with q the number of edges of the elementary polygonal tile and p 
the number of polygons meeting at each vertex. Such {p, q} tilings satisfy (p — 2)(g — 2) =4 
in the Euclidean space E2 , (p — 2)(g — 2) < 4 on 5*2 and (p — 2)(g — 2) > 4 on H2. This 
leaves the two dual triangular/hexagonal tilings {3, 6} and {6, 3} and the square tiling {4, 4} 
in E2, the five tilings {3,3}, {4,3}, {3,4}, {5,3} and {3,5}, corresponding to the platonic 
solids, in S2 and an infinity of tilings in Ho j30j . 

The fiexibility offered by the multiple tilings in H2 has for instance been used by Modes 
and Kamien[311[32] to study "isostatic" packings of hard disks in H2. Isostaticity means that 
the number of constraints coming from force and torque balance equations is exactly equal to 
the number of degrees of freedom in the system p^2] . It is a global requirement that is related 
to marginal stability in solids p^3l I104J : it has recently received renewed attention in the 
context of "jamming phenomena" present in equilibrium and driven disordered assemblies of 
particles |1U5] . Isostaticity is a topological property that does not depend on the curvature 
of space. For a two-dimensional system of N hard disks, there are 2N degrees of freedom 
and zN/2 constraints coming form the number of contacts, where z is the average number 
of contacts per disk, so that isostaticity corresponds to ^ = 4. Therefore, all allowed tilings 
{4, q} in H2 are isostatic, which broadens the scope of systems that can be used to study 
the generic features associated with isostaticity [311 1^ . 

As stressed several times in this article, only hyperbolic geometry allows one to study 
macroscopic systems in the thermodynamic limit at constant nonzero curvature. Spherical 
geometry implies a finite system which can be investigated per se, as in studies of colloidal 
systems in spherical substrates, or used as trick to converge to the thermodynamic limit in 
the Euclidean plane by decreasing the curvature. The latter approach, corresponding to the 
already mentioned spherical boundary conditions, has been taken to study dense disordered 
packings of hard particles [591 EI] as well as the nature of crystallization for particles inter- 
acting with power-law and logarithmic pair potentials [65l E^. In both cases, the rationale 
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is that spherical boundary conditions do not favor crystalhne hexagonal arrangements as 
much as periodic boundary conditions directly implemented in E2, which then provides a 
less biased finite-size approach to phenomena taking place in fiat space in the thermody- 
namic limit. On the other hand, one may be interested in studying the infiuence of geometric 
frustration on the structure and the dynamics of a system and therefore work with a con- 
stant nonzero curvature. Consider the hyperbolic plane H2 which, as already stated, is of 
infinite spatial extent. For disks of diameter a embedded in H2, the local order of the liquid 
changes as one increases the radius of curvature k,~^. From hexagonal at zero and small cur- 
vature parameter na, it becomes heptagonal at a larger na, then octogonal, etc[TT], and for 
commensurate values of the curvature parameter, the locally preferred structure can freely 
propagate to tile space; in such cases, corresponding top!T] 



K„a = 2 cosh 



2sm(^) 



(6.1) 



with n = 6, 7, 8, ..., there is no frustration. The effect of frustration can be investigated off 
these commensurate curvatures. For instance, frustrated hexagonal order can be studied for 
K(T larger than zero [hqct = 0) but significantly less than kjct ^ 1.09055 [TT]. 

Before moving on to a more extensive discussion of frustration and glass formation in 
a liquid model in the hyperbolic plane, we briefiy describe the various tools that have 
been used in practice to assess the structure of dense phases in curved space. All of them 
are actually extensions of methods and observables developed for the Euclidean case. For 
spherical particles, both "positional" (translational) and "bond-orientational" orders are 
of interest. "Positional" refers to the distribution functions introduced in Sec. HV] that 
involve correlations between the particle centers. The most easily measured or computed 
is the radial distribution function g{r), with r being the geodesic distance between two 
particle centers. "Bond-orientational" refers to the distribution functions associated with 
the (artificial) "bonds" joining two nearest-neighbor particles. There are different ways of 
defining such bonds, a commonly used one being to assign nearest neighbors through a 
Voronoi-Dirichlet construction (see below). In Euclidean space, one associates a unit vector 
to each bond and orientational order refers to the average and the correlations of local order 
parameters defined in two dimensions as [9] 
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^n(r,) = — ^exp(^n^,fc), (6.2) 

^ k=i 

where the sum is over the A^;, nearest neighbors of the particle located at r^ and 6jk is the 
angle characterizing the "bond" between atoms j and k [n = 6 for instance is character- 
istic of 6— fold hexagonal or hexatic order). Similar quantities can be introduced in higher 
dimensions |106j . 

The difficulty that one faces when dealing with non-Euclidean space is that there is 
no global existence of vector fields: vectors are defined locally (in the tangent Euclidean 
manifold) and to be combined or compared to vectors at another point in space, they 
must be "parallel transported" along the geodesic joining the points (see Appendix A and 
Refs.[7ll3B])- As a result, the bond-orientational correlation functions involving an extension 



to non-Euclidean geometry of the local order parameters in Eq. |6.2| are path-dependent in 

the presence of a nonzero curvature. A natural definition of the pair correlation functions is 

however to consider the geodesic between the two points under consideration, which leads 

to 

1 " 
^"W = N E(^"(^l^')^«(^'))r./^'Hr., -r), (6.3) 

where \Ef„(2|j) is the order parameter when parallel transported from point r^ to point r^ 
along the geodesic Fjj, rij is the geodesic distance between the two points and 5^^^ is the 
delta function that is appropriate for the non- Euclidean metric (for 5*2, see Ref.[60] and for 
H2 see Refs.|TTl|69]). Finally, the Voronoi tesselation or its dual, the Delaunay construction 
(or Dirichlet triangulation in two dimensions), is a way to uniquely define nearest neighbors, 
hence bonds, coordination number, etc, which is most useful in dense liquid and amorphous 
phases. In particular, it allows one to characterize at a microscopic level topological defects 
(especially disclinations which are orientational defects in the form of lines in three dimen- 
sions and points in two dimensions) occurring in a putative order. These constructions can 
be extended to sphericalfTJ [571 EO] and hyperbolic geometry [HI [TOl 1107] . 

B. Glassforming liquid on the hyperbolic plane 

We briefiy review here the extensive computer simulation study that we have recently 
carried out on the structure and the dynamics of the truncated Lennard- Jones liquid (see 
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(a) 



(b) 



FIG. 9: Locally preferred structure in atomic liquids: (a) Frustrated icosahedral order in 

three-dimensional Euclidean space E^: the distance I between two neighboring outer 

spheres is slightly larger than the particle diameter a which is also the distance between 

the central sphere and the 12 outer ones. |(b)| Frustrated hexagonal order in the hyperbolic 

plane H^: the nonzero curvature induces that, here too, / is slightly larger than a. 



above) embedded in H2 with curvature parameter Ka small enough that the local order in 
the liquid is hexagonal/hexatic and is therefore frustrated at long distance. The motiva- 
tion behind this study was to assess the validity of the frustration-based theory of the glass 
transition [Tool HOlj - In the latter, frustration is hypothesized to be ubiquitous in liquids. 
The salient features of the phenomenology of glassforming liquids, above all the spectac- 
ular " super- Arrhenius " increase of the viscosity and the relaxation time as one lowers the 
temperature |101j . are then attributed to the frustration-limited extension of the locally 
preferred liquid order; this extension is driven by the proximity to an "avoided ordering 
transition" that would take place in the absence of frustration |100[ 1101] . Frustrated icosa- 
hedral order in three dimensions with 5*3 as the unfrustrated spaced E] is one example that 
we have already discussed. A simpler model however is provided by frustrated hexagonal 
order in two-dimensional negatively curved space P, [TH [T2j, as illustrated in Fig. M In this 
case, "ideal ordering" takes place in flat space and this ordering transition observed at a 
temperature T* is avoided as soon as one introduces a nonzero curvature, with the curvature 
parameter na playing the role of the frustration strength. 
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FIG. 10: Slowing down of relaxation in the truncated Lennard- Jones model in H2 (MD 

simulation |68j): logarithm of the translational relaxation time r versus T* /T for p ^ 0.85 

and for various curvature parameters na (T* is the ordering temperature in E2). The 

dotted line is the Arrhenius T-dependence roughly observed at temperatures above T*. 

When KO" > 0, the system remains liquid below T* until it forms a glass. The deviation 

from Arrhenius behavior increases when curvature (hence frustration) decreases. 



Glass formation which is preempted by ordering at T* in the Euclidean plane becomes 
possible in the hyperbolic plane and the liquid phase can be kept in equilibrium at temper- 



atures below T* . This is shown in Fig. 10 where we plot the logarithm of the relaxation 
time versus 1/T for several curvatures, as obtained from Molecular Dynamics simulation. 
(Details on the Molecular Dynamics simulation technique in hyperbolic geometry and on 
the generalization of the time-dependent correlation functions are given in Refs. [0^1470] .) 
One observes that curvature plays virtually no role for temperatures above T* and that 
super-Arrhenius behavior, i.e. a deviation from simple Arrhenius T— dependence, becomes 
significant around T* and is more pronounced as one decreases the curvature parameter na. 
These observations support the theoretical predictions that the avoided transition controls 
the slowing down of the relaxation and that the "fragility" of a glassformer, which quantifies 
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how much it deviates from simple Arrhenius behavior, decreases as one increases the frus- 
tration, i.e here the curvature. Note that this trend should also apply to liquids of colloidal 
particles on a spherical substrate: at least for small enough curvature such that the local 
order remains hexagonal, decreasing the curvature should lead to a stronger slowing down 
of the dynamics. 

The effect of the curvature can also be seen on the structure of the liquid. This is 



illustrated in Fig. 11 where the radial distribution function g{r) of the Lennard- Jones model 
in the Euclidean plane is compared to that in H2 with na = 0.1. The (quasi) long-range 
order present below T* in fiat (unfrustrated) space is no longer observed in curved space, 
thereby confirming that the one-component system remains a disordered liquid in the latter. 
On general grounds, one expects that curvature-generated frustration induces three dif- 
ferent regimes in a liquid as temperature decreases. In a first regime, above T*, the structure 
and the dynamics only involve local properties; as, locally, curved space looks fiat, there is 
no significant influence of curvature (at least for small curvatures that accommodate the 
same local order). A second regime near and below T* is controlled by the proximity of the 
avoided transition: the locally preferred structure extends in space and the associated corre- 
lation length grows with decreasing temperature (as can be checked by studying the bond- 
orientational order correlation function G^lr) = g%{T)/g{r), see Ref. [69]). A final regime 
is reached when the correlation length associated with growing (bond-orientational) order 
saturates due to frustration: the spatial extent of the frustrated order cannot grow beyond 
the radius of curvature and the final regime is dominated by the presence of an irreducible 
density of topological defects. More will be said below about this regime. This three-regime 
scenario has been observed in computer simulations of liquids in both negative [69l [70] and 
positive [60] curvature manifolds. 

VII. GROUND-STATE PROPERTIES, ORDER AND DEFECTS 

So far, we have mostly dealt with liquid and fluid phases in curved space. We briefly 
discuss now low-temperature (or high-density) phases and ground-state properties. In the 
preceding section, we have stressed that a nonzero curvature can either induce long-range 
order |1 08] . allowing tiling of a homogeneous curved space by regular polytopes for specific 
values of the curvature, or, on the contrary, frustrate the long-range order present in Eu- 

49 



5 

4 

_^ 3 




' y 






' 






:■ 








:■ 


DO 








L 


2 








r 









': 


1 


















i 



5 10 

r/a 

(a) 



15 




r/a 



10 12 14 



(b) 



FIG. 11: Radial distribution function of the truncated Lennard- Jones model in 7^2 • The 
(blue) curves correspond to temperatures above the flat-space ordering transition 



(T > T*). (a) Euclidean plane {na = 0): when T <T* (red curves) g{r) displays many 

peaks corresponding to the existence of a quasi-long range positional order. |(b) Hyperbolic 

plane with na = 0.1: when T <T* (red curves), there is no signature of quasi-long range 

order and the system remains in a liquid phase. 
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clidean space. We focus on the latter situation and we moreover restrict the discussion to 
the two-dimensional case with spherical particles, in which the hexagonal/triangular order 
that forms the ground-state in flat space is frustrated by the introduction of curvature. The 
topic of order, curvature and defects in two dimensions has recently been nicely reviewed 
in a quite exhaustive article by Bowick and GiomipT]. Therefore, we only intend to give a 
brief survey of this aspect, for the sake of completeness of this article. 

A first insight into frustrated hexagonal order in non-Euclidean geometries is provided 
by topological considerations. Consider the Delaunay triangulation of a dense assembly 



of particles on a two-dimensional manifold of genus g (see Sec. VI) and, if open, with h 
boundaries. Its Euler characteristic is then x = 2(1 — g) and the Euler-Poincare theorem 
states that any "tesselation" (i.e. tiling) of the manifold satisfies 

V-E + F = x, (7.1) 

where V, E and F are the numbers of vertices, edges and faces in the tesselation. In the 
case of the Delaunay tesselation by triangles, this leads to 

N 

-^(6-^) = X, (7.2) 

where iV = V^ is the number of particle centers (i.e., of vertices) and z is the average 
coordination number of the particles. If one defines the topological charge g^ of a particle 



with coordination number Zj as g^ = 6 — Zj, one can rewrite Eq. |7.2|as 



N 

J]g, = 6x = 12(l-^)-6/i. (7.3) 

For the Euclidean plane, g = 1 and /i = 0, so that x = ^"^^ z = 6: perfect hexagonal tiling 
is possible. For the sphere, 5*2, g = 2 and h = 0, so that the average coordination number 



must be strictly less than 6. The minimal way to satisfy Eq. |7.3| is then to have 12 particles 
with configuration number number 5 in an otherwise 6— fold coordinated configuration. Such 
particles appear as (point) topological defects in hexagonal order. They represent positively 
charged (g = +1) disclinations. For the hyperbolic plane H2, one may first take a detour 
via periodic boundary conditions (see Appendix B). When imposing a periodic boundary 
condition, the primitive cell containing the system corresponds to a (compact) quotient space 
of genus g > 2, hence with x = —2{g — 1). In consequence, the total topological charge 
in the primitive cell has to be equal to —I2{g — 1): there must be an excess of negatively 
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charged disclinations in this case (e.g. of dischnations with q = —1 corresponding to 7-fold 
coordinated particles). By using geometric input in the form of the Gauss-Bonnet theorem 
that relates curvature K = —k? and characteristic x, 

JdSK = 27ix, (7.4) 



where we recall that dS = ^/\g{x)\d'^x (see Eq. (4.2)) and S is the primitive cell (or funda- 
mental polygon) of the periodic boundary condition, one finds that the area of the primitive 
cell is given hj A = —27txk,^'^- As a result, Eq. (7.3) becomes 



lf,. = _(!).^ (7.5) 



i=l 



which means that the density of topological charge in H2 is only controlled by the curvature, 
irrespective of the choice of periodic boundary condition, thereby providing a well defined 



thermodynamic limit for a "bulk" property (see also the discussion in Sees. Ill and IV). 

In curved two-dimensional space, hexagonal order must come with an irreductible number 
(or density) of topological defects even in the ground state. Actually, constructing the 
ground state of a system of particles on a two-dimensional manifold is highly nontrivial. 
The question goes back to Thomson who considered the ground state of repulsive charged 
particles on a sphere p^9] . (Another famous example is the explanation by Caspar and 
Klug |llG] of the icosahedral symmetry of spherical virus capsids.) Brute- force numerical 
approaches are plagued by the presence of multiple low-energy minima, which are induced 
by geometric frustration, and finding the global minimum for a large number of particles is 
extremely difficult. To get around this problem, coarse-grained approaches have been devised, 
most prominently a continuum elastic theory developed by Nelson, Bowick, Travesset and 
their coworkers pTf 111111116 ]. In the latter, one directly deals with the defect degrees of 
freedom, and all the microscopic information about particle interactions is embedded in 
effective elastic constants and defect core energies j57]. (Note that an alternative route to 
the continuum elastic theory has recently been proposed on the basis of a coarse-grained 
density functional theory: see Ref. |117j .) 

In the Bowick- Nelson- Travesset approach |112j . one focuses on the density of disclina- 
tions, which are the elementary defects from which other defects such as dislocations can be 
built, and treats the 6-fold coordinated particles through continuum elastic theory. At low 

52 



temperature, the elastic free energy of an arbitrary disclination density 



( \ ^ 

s X 



where iV^ is the total number of disclinations, is given by [271 1112] 

(7.7) 
where F is the Young modulus of the hexagonal crystal in flat space, Ecore,i is a disclination- 
core free energy renormalized by thermal fluctuations (112] . i^(x) is the Gaussian curvature 
at point X, and 6*2 (x,y) is the Green function of the bi-Laplacian (or biharmonic operator) 
on the manifold S, i.e. satisfying 

A2G2(x) = 5(2)(x), (7.8) 

with A the Laplace-Beltrami operator. The first term of the free-energy functional is min- 
imized by having s(x) = -ft'(x) at each point, which means that the distribution density 
exactly cancels the effect of the Gaussian curvature. If one considers homogeneous mani- 
folds of constant Gaussian curvature, K{x.) = K, such as 5*2 and H2, the former condition 
cannot be everywhere satisfied as the disclinations are discrete objects. Topology however 
requires a generalization of the electroneutrality constraint, namely 

J d'x./W)\isi^)-K)=Q, (7.9) 



which by inserting Eq. (7.6) and the Gauss-Bonnet theorem, Eq. (7.4), is equivalent to 
Eq. (7.6). Therefore, on average, the disclination charges screen the Gaussian curvature. 

As the core energy of the dislocations increases with the magnitude of their topological 
charge, the second term of the free-energy functional is minimized by having the smallest 
irreducible number of disclinations with elementary topological charges, e.g. 12 disclinations 
of charge g = +1 on 6*2 and a density of disclinations of charges g = — 1 equal to -k^ on 7^2- 

For small enough curvature (which for 5*2 means a large enough number of particles if the 
particle size a is kept fixed), adding extra disclinations in the system on top of the irreducible 
number may lower the (free) energy by screening more efficiently the curvature [27 1 11121 1118] . 



The total charge of these extra disclinations is zero (so that Eq. (7.3) remains satisfied) and 



the most economical way for them to organize is by forming dislocations that consist of 
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(a) 



(b) 



FIG. 12: Illustration of grain boundary scars in a hexagonal background in 
two-dimensional curved space, (a) Model grain boundaries in 5*2 obtained by minimization 

of Eq. (7.7) [IH]. (b) Low-temperature atomic configurations of the truncated 
Lennard- Jones liquid in the Poincare disk representation of if ^[511 (see Appendix A for the 

definition of Poincare disk) . 



"dipoles" made by a positive and a negative disclination. The surprising outcome of the 
theoretical studies on order on the sphere 5*2 is that these dislocations form strings that 
radiate from each of the irreductible disclinations and terminate in the hexagonal crystalline 
background at a finite distance [23 11121 I118J . These structures which are found found on 
the hyperbolic plane [7D], but are forbidden in Euclidean space, have been dubbed "grain 



boundary scars "[IHl I112J . We illustrate these structures on 5*2 [18] and if^2[Zn] in Fig. 12 
The theory predicts the number of dislocations per irreductible disclination as well as the 
cut-off distance Vc before which the strings of dislocations terminate. For instance, the latter 
is given in 5*2 by |112] 



and in if 2 by [70 



R arccos 



K ^arccosh 



0.59i?, 



0.57i?, 



(7.10) 



(7.11) 



regardless of the microscopic details of the system of particles. These predictions concerning 
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grain boundary scars in 5*2 and H2 are well supported by experiment alp ^ [1 19] and computer 
simulation results [23, EOl 1117] . The elastic theory based on defects has been generalized to 
long range interactions between particles on a sphere, to other two-dimensional manifolds, 
with varying curvature and possibly open, and to other type of order: this is discussed in 
detail in the review of Bowick and Giomi[27]. where all relevant references can also be found. 
In addition, the dynamics of the defects, essentially the diffusion of the dislocations, has been 
studied theoretically and experimentally on the sphere [T71 [271 1113] as well as theoretically 
and by computer simulation in the hyperbolic plane [70]. 

VIII. CONCLUSION 

In this article, we have reviewed the work on the structure and the dynamics of fluids, 
liquids, and more generally dense phases, in curved space. We have stressed that the moti- 
vation for such studies is twofold. On the one hand, curving space provides an additional 
control parameter, curvature, to shed light on the behavior of systems of interest in the fiat, 
Euclidean space. Examples are provided by the use of spherical and hyperspherical bound- 
ary conditions as well as by investigations of the glass transition in the context of geometric 
frustration. On the other hand, there are many physical situations in physical chemistry, 
soft condensed matter and material science in which a two-dimensional curved surface is 
coated by a layer of particles that can move and equilibrate on the curved substrate. We 
have focused on cases where the geometry of the substrate is frozen, and most theoretical 
developments have been presented for the case of homogeneous two-dimensional manifolds 
of constant Gaussian curvature, the sphere 5*2 and the hyperbolic plane H2. 

Non-Euclidean geometries bring in subtleties in the statistical mechanics of systems of 
particles when compared to the standard Euclidean case. We have discussed the main ones, 
concerning the thermodynamic limit and the effect of the boundaries, the definition of the 
pressure and the relations involving correlation functions, the constraints on the form of the 
interaction potentials, the properties of particle diffusion, or the nature and the organization 
of the defects in dense and quasi-ordered phases. Additional results can be summarized as 
follows: curvature has a negligible or weak effect on the fluid behavior at high temperature 
and/or low density when the local structure or dynamics are probed. On the other hand, 
a nonzero curvature has strong consequences in dense phases where it can either induce 
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or frustrate ordering. In the latter case, it slows down the relaxation and allows glass 
formation, as well as imposes in an irreducible number of topological defects in ground-state 
and low-temperature configurations. Curvature has also a distinct influence when long- 
range interactions, as in Coulombic systems, or long-range correlations, as in the vicinity of 
a critical point, are expected: there, the radius of curvature acts as a cutoff that imposes 
exponential decay at long distances. 

Needless to say that progress should be made to describe liquids and fluids on substrates 
with varying curvature or even with fluctuating curvature as encountered in membranes. 
Above all, one can hope that in the near future more experiments will be performed on 
curved substrates to provide systematic information on the phase behavior and the dynamics 
of particle systems for a significant range of temperature or density and for several curvatures 
or geometries. 

Appendix A: A recap on Riemannian manifolds 

To describe liquids on curved surfaces, let us introduce some elements of differential 

geometry [3 133 EH]- A rf— dimensional Riemannian manifold is defined by a set of coordinates 

X = (xi,a;2, ■■■,Xd) and a metric tensor Qiji^s). The length element ds is defined by 

d 
ds^ = y^gij{'x.)dxidxj, (Al) 

and the "volume" element is given by 

d 

dS = ^/W)\l[dx,, (A2) 

where |fi'(x)| denotes the absolute value of the determinant of the metric tensor. 

Diagonalizing the curvature tensor provides the principal directions of the manifold 
(eigenvectors) at point x and the eigenvalues are the radii of curvature, Ri,i = l,...,d. 
The mean curvature is defined as the algebraic average of the curvatures. 

In the case where the manifold is two-dimensional, one defines the Gaussian curvature as 
the product of the two curvatures, 

KU) = — , } , , . (A4) 

^ ^ i?i(x)i?2(x) ^ ' 
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*:.„(A-)^^^E^™^ (A5) 



igradif)y = J29''^^ (A6) 



If X is a vector field (defined at each point x in the tangent Euchdean manifold), the 
divergence is given by 

and the gradient of the scalar function / 

d 

where g^^i'x) is the inverse tensor of gij{^), i.e. Ylij=i9i39^^ — ^i with 6^ the Kronecker 
symbol. 

The Laplace-Beltrami operator A acts on a function /(x) as A/ = div grad{f), i.e., 

We mostly consider in this article homogeneous manifolds of constant nonzero curvature. 
In two dimensions, there are two associated geometries, the spherical and the hyperbolic 
ones. For the sphere 5*2, the metric tensor is diagonal in the angular coordinates 9 (colati- 
tude) and (p (longitude) and one has 

gee = R\ g^^ = R^sm{9f. (A8) 

The squared length element ds'^ is then equal to 

ds^ = R\de^ + sm{efd(t)''), (A9) 

the area element is given by 

dS = R^ sm{e)d(j)de, (AlO) 

and the two radii of curvature are equal and constant with Ri = R2 = R and the Gaussian 
curvature K = R~^. Finally, the Laplace-Beltrami operator is 

For the hyperbolic plane H2 (also called "pseudosphere " or "Bolyai-Lobachevski plane ") [291 
130] , one can use the polar coordinates r and 0. The metric tensor is then diagonal with 
grr = 1 and g^^p = K^^sinh^KrY, which gives a squared length element 

.2_ > 2 , /sinh(fi:r)\ ^ 



ds' = dr'+[ ^^ ) d(P' (A12) 
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and an area element 



dS = ''^^^^'\ rd<p. (A13) 



The two radii of curvature are of opposite signs, Ri = —R2 = n ^, so that K = —k^. In 
addition, the Laplace-Beltrami operator is given by 

A = ^ /^_^sinh(Kr)— ^ + ( '^ \ — ■ (A14) 

sinh(fi;r) \dr dr J \smh(Kr) J dcj)^ 

The hyperbohc plane cannot be embedded in three-dimensional Euclidean space (contrary 
to S2) and "models", i.e. projections must be used for its visualization. A convenient one is 
the Poincare disk model which projects the whole hyperbolic plane H2 onto a unit disk. The 
projection is conformal (angles are conserved) but not isometric (distances are deformed and 



shrinks as one reaches the disk boundary). This representation is used in Fig. 12 d. If x and 
y are the Cartesian coordinates of a point on the unit disk, the relation to the above polar 
coordinate is as follows: 



V^^T^ = tanh (^y j , (A15) 

- = tan(0), (A16) 

X 

and the squared length element ds^ is given by 

d,= = K-^M±M_, (A17) 

(1 - [;x^ + v-'-)f 

From the above formulae, one can compute for instance the geodesic distance ri2 between 
two points ri and Y2- In 6*2, one finds 

cos(6'i2) = cos(6'i) cos(6'2) - sin(6'i) sin(6'2)cos(0i - ^2) (A18) 

where (6'i,0i) and (6'2,02) are the coordinates of the two points and R612 the geodesic 
distance between these points. On the other hand in H2., one has 

cosh(Kri2) = cosh(Kri) cosh(/T;r2) — sm\i{nri) smh{nr2) cos(0i — 02), (^19) 

where (ri, <^i) and (r2, (^2) are the (polar) coordinates of the two points. Note the symmetry 



between Eqs. (A18) and (A19) in the exchange of i? -H- m ^ . The same symmetry takes 



place in relating spherical and hyperbolic trigonometries. So for instance, from the standard 
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result on a sphere, one finds the following trigonometric relations for the a general hyperbolic 
triangle with sides a, b, and c and opposite angle a, (5 and 7(30]: 

sinh(/T;a) smh{Kb) sinh(/T;c) 
sin(a) sin(/3) sin(7) 

cosh(/T;c) = cosh(Ka) cosh(/T;6) — sinh(Ka) s\.r\h{td)) cos(7), (A21) 

cosh(.c) = ^os(a) cos(/3) + cos(7) _ ^2) 

sin(a) sin(p) 

Note that the first two relations have Euclidean counterparts (obtained by letting k — > 0), 
but not the third one which is specific to nonzero curvatures. 

Appendix B: Periodic boundary conditions on the hyperbolic plane 

Generically, implementing periodic boundary conditions consists in choosing a primitive 
cell containing the physical system such that it can be infinitely replicated to tile the whole 
space. So, prior to building periodic boundary conditions, one needs to know the allowed 
tilings of the space under consideration. Here, for simplicity, we will limit ourself to regular 
cells and so to regular tilings. 

On the hyperbolic plane if 2, an infinity of regular tilings {p-,q\ are allowed if p (the 
number of edges of the primitive cell) and q (the number of cells meeting at each vertex of 
the tiling) verify the following condition 

(p-2)(g-2)>4. (Bl) 

This, therefore, opens the possibility to have an infinite number of possible periodic boundary 
conditions. 

To ensure smoothness and consistency, the edges of the primitive cell of any periodic 
boundary condition have to be paired in a specific way: leaving the cell through one edge 
implies to come back by another edge, a process which should be equivalent to exploring 
the whole tiling of the plane. Constraints thus arise on how cell edges are paired together 
in addition to those on the shape of the cell. We give in the following a rapid overview of 
how to classify and construct periodic boundary conditions in H^ by describing cell shapes 
and edge pairings. 

First, a fundamental polygon (primitive cell with properly paired edges) encodes all the 
needed information to build periodic boundary conditions and to replicate the system in 
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FIG. 13: Simplest fundamental polygon on the hyperbolic plane. It is associated to the 
{8, 8} tiling. The arrows indicate the way edges are paired. 





(a) 



(b) 



FIG. 14: Schematic representation of the "compactification" of the fundamental polygon 



shown in Fig. 13 The paired edges are glued together: 1 with 3 and 5 with 7 in (a); 2 with 
4 and 6 with 8 in (b) . The final compact manifold is a two-hole torus represented in 



Fie. 15 



the entire space. The simplest (and smallest) fundamental polygon in H^ is an octagon, 
corresponding to an {8, 8} tiling, with the edge pairing shown in Fig. [13] By gluing the 



paired edges together (see Fig. 14 for of visualization of intermediate states) it is possible 
to represent the fundamental polygon as a compact manifold, also corresponding to the 
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FIG. 15: Compact manifold and graph obtained by gluing the edges of the fundamental 

polygons of Fig. 



13 



"quotient space" (see Fig. 15). The above octagonal periodic boundary condition leads to a 
2-hole torus, whose genus (number of holes, here g = 2 ) fixes the area of the fundamental 



polygon through the Gauss-Bonnet theorem (see Eq. (7.4): A = Ann {g — 1)). As the 



genus is an integer, compact manifolds embedded in the hyperbolic plane can only lead to a 
discrete set of areas. In more technical words, an homothety cannot be applied to hyperbolic 
manifolds without changing the curvature. Here, the octagonal fundamental polygon cannot 
be scaled at constant curvature to allow one to study a bigger system for instance. Therefore, 
to change the area of the fundamental polygon in H^, one as to vary its genus, which in 
turn implies to change its symmetry, more particularly its number of edges. 

To classify all the possibilities (for regular polygons), one can use the properties of the 
graph formed by glued edges of the fundamental polygon and embedded in the associated 



(yf-hole torus (see Fig. 15 for an example on the 2-hole torus). The constraints on such 
graphs to obtain relevant periodic boundary conditions are detailed in [M] which allows 
one to classify and build all possible regular periodic boundary conditions in the hyperbolic 
plane |120] . The classification exhibits "families" that comprise graphs with a given number 
of vertices and a given pairing pattern but different values of the genus g and that share 
similarities when varying g. In this framework, the octagonal fundamental polygon shown 



in Fig. 13 can be seen as the direct generalization of the square periodic boundary condition 



encountered in the Euclidean plane (both are in the same "family"). 
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